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ABSTRACT 

As a part of the Herschel key programme PRISMAS, we have used the Herschel-HWl instrument to observe interstellar nitrogen 
hydrides along the sight-lines towards eight high-mass star-forming regions in order to elucidate the production pathways leading 
to nitrogen-bearing species in diffuse gas. Here, we report observations towards W49N of the NH N = 1-0, 7 = 2-1, and 7=1-0, 
ortho-NH 2 Nk^jCc J = 1 1,1 3/2 - Oo.o 1/2, orfho-NH 3 J K = l -0 and 2 - lo, para-NH 3 J K = 2i — lj transitions, and unsuccessful searches 
for NH + . All detections show absorption by foreground material over a wide range of velocities, as well as absorption associated 
directly with the hot-core source itself. As in the previously published observations towards G10. 6-0.4, the NH, NH 2 and NH 3 
spectra towards W49N show strikingly similar and non-saturated absorption features. We decompose the absorption of the foreground 
material towards W49N into different velocity components in order to investigate whether the relative abundances vary among the 
velocity components, and, in addition, we re-analyse the absorption lines towards G10. 6-0.4 in the same manner. Abundances, with 
respect to molecular hydrogen, in each velocity component are estimated using CH, which is found to correlate with H 2 in the solar 
neighbourhood diffuse gas. The analysis points to a co-existence of the nitrogen hydrides in diffuse or translucent interstellar gas 
with a high molecular fraction. Towards both sources, we find that NH is always at least as abundant as both o-NH 2 and o-NH 3 , 
in sharp contrast to previous results for dark clouds. We find relatively constant A'(NH)/A'(o-NH 3 ) and _/V(o-NH 2 )//V(o-NH 3 ) ratios 
with mean values of 3.2 and 1.9 towards W49N, and 5.4 and 2.2 towards G10.6-0.4, respectively. The mean abundance of o-NH 3 
is ~2xl0~ 9 towards both sources. The nitrogen hydrides also show linear correlations with CN and HNC towards both sources, and 
looser correlations with CH. The upper limits on the NH + abundance indicate column densities <2- 14% of A'(NH), which is in 
contrast to the behaviour of the abundances of CH + and OH + relative to the values determined for the corresponding neutrals CH and 
OH. Surprisingly low values of the ammonia ortho-to-para ratio are found in both sources, w 0.5-0.7+0.1, in the strongest absorption 
components. This result cannot be explained by current models as we had expected to find a value of unity or higher. 

Key words. ISM: abundances - ISM: molecules - Sub-millimetre: ISM - Molecular processes - Line: formation - Astrochemistry 



1. Introduction 

Nitrogen is among the six most abundant elements in the uni- 
verse and, despite its fundamental role in the chemistry of 
molecules connected with life, the chemical network of nitro- 
gen in the interstellar medium is still poorly understood due to 
severe difficulties to observe key molecules from the ground. 



Send offprint requests to: carina.persson@chalmers.se 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



Today about 55 molecules containing nitrogen have been dis- 
covered in interstellar space and a few more in circumstellar en- 
velopes 1 . The major reservoir of nitrogen is believed to be in 
the form of atomic N and molecular N2. The latter is extremely 
difficult to observe since it has no permanent dipole moment 
(Sandford et al. 2001). In dense molecular clouds N2H + is in- 
stead often used as a tracer of N2, which remained undetected 
until Knauth et al. (2004) reported far-ultraviolet observations 
towards HD 124314. The total abundance of nitrogen therefore 



1 www.astrochymist.org, www.astro.uni-koeln.de/cdms/molecules 
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still relies, to a high degree, on chemical modelling and observa- 
tions of nitrogen-bearing species other than N2. 

In order to constrain the nitrogen formation pathways, obser- 
vations of nitrogen hydrides are crucial since they are at the root 
of the nitrogen chemical network, appearing in its first steps in 
chains of reactions that lead to other more complex species. The 
abundances of these species are thus key diagnostics for the ni- 
trogen chemistry. The nitrogen hydrides are, however, also prob- 
lematic to observe since their ground state rotational transitions 
lie at sub-mm wavelengths and are thus very difficult, or impos- 
sible, to reach from the ground. Key species, such as imidogen 
(NH) and amidogen (NH2), have therefore previously not been 
widely observed, and there is still no detection of the NH + radi- 
cal. 

Although very few observations exist of NH and NH2 in in- 
terstellar space, they are well known in comets (e.g. Swings et al. 
1941; Meier et al. 1998; Feldman et al. 1993), and have been 
observed in stellar photospheres (e.g. Schmitt 1969; Farmer & 
Norton 1989) via their electronic, vibration-rotation, and high 
rotational transitions. The first detection of interstellar NH was 
made by Meyer & Roth (1991) by optical absorption spec- 
troscopy. Subsequent observations by Crawford & Williams 
(1997) and Weselak et al. (2009) have yielded several lines of 
sight in diffuse and translucent gas where column densities of 
NH, CH, CN, and H2 have been directly measured. The average 
value of the column density ratio in these diffuse and translucent 
sight-lines is /V(NH)//V(H 2 ) = 3xl0~ 9 . 

Interstellar NH2 was first detected by van Dishoeck 
et al. (1993) in absorption towards SgrB2 in three fine- 
structure components of the para-symmetry form of NH2, the 
Nk„,k c = 1 1,0 _ lo,i transition, with partially resolved hyperfine 
structure at frequencies 461 to 469 GHz. The Infrared Space 
Observatory (ISO) was later used to observe unresolved absorp- 
tion lines of both ortho and para symmetry forms of NH2, as 
well as NH, towards this source through the use of the long- 
wavelength spectrometer (Cernicharo et al. 2000; Goicoechea 
et al. 2004; Polehampton et al. 2007). 

In contrast to NH and NH2, ammonia (NH3) has been ex- 
tensively observed for more than 40 years and was in fact the 
first polyatomic molecule to be identified in interstellar space 
(Cheung et al. 1968) by means of its J = K inversion transi- 
tions at cm wavelengths (K is the quantum number of the pro- 
jection of total angular momentum J on the molecule's sym- 
metry axis). The ammonia molecule has, however, two sym- 
metries, like NH2, which arise due to the possible orientations 
of the hydrogen spins and behave like two distinct species: 
ortho-NH 3 (all H spins parallel, K-3n where n is an integer 
> 0) and para-NH 3 (not all H spins parallel, K + 3n). Important 
information about the ammonia formation pathways could be 
inferred from observations of both symmetries. Unfortunately, 
only the para form has relatively low-excitation transitions ac- 
cessible from ground since the K-0 ladder of energy levels 
has no inversion splitting and the Jk — 33 inversion transitions' 
lower energy level is 122K above ground. Ortho inversion lines 
(K = 3n, n > 1) can thus only be observed in relatively warm 
molecular gas. Only a few ammonia observations of the cold dif- 
fuse interstellar gas exist, using para inversion lines (Nash 1990; 
Tieftrunk et al. 1994; Liszt et al. 2006) which leaves the am- 
monia formation mechanism poorly constrained in diffuse gas. 
The (0,0) ammonia ground state, with ortho symmetry, can only 
be studied by sub-mm rotational transitions. Observations of 
the fundamental rotational transition of ortho-NH3 J K = lo-Oo 
at 572 GHz, which has similar upper state energy as the in- 
version lines, but several orders of magnitudes higher critical 



densities, thus requires space telescopes. The Kuiper Airborne 
Observatory (Keene et al. 1983) performed the first observa- 
tions of the J K = lo-Oo transition using heterodyne receivers, 
and later on the Odin satellite continued such observations in 
a number of environments, for instance in photo-dissociation 
and star-forming regions (e.g. Larsson et al. 2003; Persson et al. 
2009), circumstellar envelopes (Hasegawa et al. 2006), diffuse 
clouds in absorption towards SgrB2 (Wirstrom et al. 2010), and 
in comets (Biver et al. 2007). 

With the launch of the Herschel Space Observatory (Pilbratt 
et al. 2010) in May 2009, unique opportunities to perform obser- 
vations of transitions between 157 and 625 /im (0.48- 1.9 THz) 
became feasible with the Heterodyne Instrument for the Far- 
Infrared (HIFI; de Graauw et al. 2010) owing to its very high 
sensitivity and spectral resolution. This allowed, for the first 
time, searches for spectrally resolved, rotational transitions in- 
volving the ground states of NH + , NH, NH 2 , and NH 3 with 
the same instrument. Several observations have already been 
reported. For instance, using Herschel-HIFl, Bacmann et al. 
(2010) found very high column densities of NH and ND in 
the cold envelope of the class protostar IRAS 16293-2422 
(2xl0 14 cirT 2 and ~1.3xl0 14 cirT 2 , respectively). Hily-Blant 
et al. (2010) found NH:NH 2 :NH 3 abundance ratios of -5:1:300 
towards the same source. 

The PRISMAS 2 key programme (PRobing InterStellar 
Molecules with Absorption line Studies) is targeting ab- 
sorption lines in the line-of-sight towards eight bright 
sub-millimetre-wave continuum sources using Herschel-HIFl: 
G10.6-0.4 (W31C), W49N, W51, G34.3+0.1, DR21(OH), 
SgrA (+ 50 km s" 1 cloud), G005.9-0.4 (W28A) and W33A. 
High-resolution absorption line spectroscopy is generally a very 
sensitive and model-independent method for measuring column 
densities of interstellar molecules, and a powerful tool to probe 
the diffuse interstellar gas clouds with no or little excitation. 

The first results and analysis of absorption lines of nitrogen 
hydrides along the sight-line towards the massive star-forming 
region G10. 6-0.4 (W31C) have already been presented in 
Persson et al. (2010, hereafter paper I). Similar abundances with 
respect to the total amount of hydrogen /V H = 2 /V(H 2 )+/V(H), 
were found for all three species: approximately 6xl0~ 9 , 3xl0~ 9 , 
and 3 x 10 9 for NH, NH2, and NH3, respectively. They were 
estimated across the whole line-of-sight and using the high 
temperature ortho-to-para limits of three and one for NH2 and 
NH3, respectively. NH + was not detected at a lcr rms level 
of 74 mK with a resolution of 1.1 MHz. The abundance pat- 
terns that we see in diffuse molecular gas are thus clearly very 
different from those in IRAS 16293-2422 and Sgr B2, where 
NH:NH 2 :NH 3 -1:10:100 and the fractional abundance of NH 
is a few times 10~ 9 . The Sgr B2 results may, however, not be 
representative of cold dark clouds since this source is very com- 
plex and atypical. 

The unexpectedly high NH abundance has been difficult to 
explain with chemical models and both the NH and NH2 pro- 
duction by purely gas-phase processes is inhibited by a lack of 
a sufficient source of N + . The models fail to simultaneously 
predict the absolute and relative abundances of the nitrogen 
hydrides. Typical steady state dark cloud models (n = 1x10 - 
5xl0 4 cm 4 , r = 10-40 K, A v > 10), predict an NH 2 abun- 
dance of (1-I0)xl0~ 8 , an NH abundance 10 times lower, and an 
NH2/NH3 ratio of 0.3 -1.5 for a wide range of assumptions (e.g. 
Langer & Graedel 1989; Millar et al. 1991), but these are not 
directly applicable to diffuse molecular gas. Examples of chem- 



2 http://astro.ens.fr/7PRISMAS 
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Table 1: Observed transitions towards W49N and G10.6-0.4. 
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416 
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AT, = 1 2 <- 0i 
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3.9 


0.018 


0.018 


0-NH2* 


952.578 


N Ka ,K c J= l u 3/2 ^ ,ol/2 


3b 
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92 


68 


2.6 


3.6 


0.018 


0.017 


o-NH 3 ' 


572.498 


Jk = lo «- Oo 


lb 


87 


1024 


94 


0.61 


0.93 


0.013 


0.025 




1214.859 


Jk = 2 q <— 1 


5a 


1024 


293 


196 


3.4 


5.2 


0.032 


0.025 


p-NH 3 


1 215.246 


Jk = 2i «- li 


5a 


1024 


293 


196 


3.4 


5.2 


0.032 


0.025 


NH + 


1012.540 


2 n 1/2 iV=l^l 7 = 3/2^1/2 


4a 


327 


171 


91 


3.0 


4.4 


0.013 


0.013 



Notes. (a) The frequencies refer to the strongest hyperfine structure components. {b) HIFI consists of 7 different mixer bands and two double 
sideband spectrometers. All transitions were observed in the upper sideband except NH + . <c) System temperature. (d) The averaged on-source 
integration time for each transition. (e) The single sideband (SSB) continuum intensity (the observed double sideband (DSB) continuum divided 
by two). ^ The rms noise, at a resolution of 1.1 MHz, divided by Tc- <s) The quantum numbers for the rotational transition Nj are F\=1 H + J and 
F = I N + F N (Klaus et al. 1997). <A) The quantum numbers for the rotational transition N KatKc J are F t =1 N + J and F = I H + Fi (Miiller et al. 1999). 
(,) The rotational energy is given by the two principal quantum numbers (J, K), corresponding to the total angular momentum and its projection 
along the molecular axis (e.g. Ho & Townes 1983; Rist et al. 1993). 



ical models for diffuse cloud conditions are found in Fig. A.l 
and A. 2 in paper I. These models were also unable to explain the 
observed abundances and ratios. Processes on dust grains have 
previously been proposed as a way to increase the NH produc- 
tion (Meyer & Roth 1991; Wagenblast et al. 1993). Such models, 
however, often predict up to 1 000 times more NH3 than NH2 
(Hasegawa & Herbst 1993, paper I). The importance of grain 
surface chemistry in diffuse clouds is also not clear since wa- 
ter ice mantles have not been detected in diffuse gas, and strong 
(UV) radiation fields counteract molecular formation on grains. 
Grain surface production of NH is therefore still debated, and 
would, if true, change our understanding of surface chemistry in 
diffuse gas. On the other hand, if grains indeed were unimpor- 
tant in diffuse gas nitrogen chemistry it would imply that either 
key gas-phase reactions must have been overlooked, or that the 
uncertainty of some reaction rates could make a difference. Both 
additional high-quality observations and chemical modelling are 
needed to solve this problem. 

In this paper, we present new observations and analyses 
of absorption lines in the line-of-sight towards the high-mass 
star-forming region W49N, and, we also re-analyse the absorp- 
tion towards G10. 6-0.4 in more detail. W49N is one of the 
most luminous high-mass star-forming regions in the Galaxy 
(~1O 7 L ) with a core that contains more than a dozen ultra- 
compact H II regions (Dreher et al. 1984; Dickel & Goss 1990; 
De Pree et al. 2000). It is located on the far side of the Galaxy 
at a distance of 11.4kpc with Galactic coordinates / = 43.17° 
and b = 0.012°, in one of the most massive giant molecular 
clouds (~ 10 6 M ) in the Milky Way. The source velocity is about 
+8kms _1 and the foreground gas along the line-of-sight is de- 
tected at vlsr ~ 30 - 75kms _1 , revealing gas at two locations 
in the near and far side of the Sagittarius spiral arm around 
40 and 60kms _1 (Dame & Thaddeus 1985). The ultra-compact 
HII region G 10.6-0.4 in the star-forming W31 complex is 
an extremely luminous sub-millimetre and infrared continuum 
source. The source is located within the so-called 30km s _1 arm 
at a kinematic distance of 4.8 kpc (Fish et al. 2003). The gas 
associated directly with G10. 6-0.4 is detected at a systemic 
source velocity of v L sr ~ -1 km s _1 , determined from OH maser 
emission observations, while the foreground gas is detected at 
vlsr ~ 10-55 kms -1 . 



Section 2 summarises the observations and data reduction, 
and the results from our Herschel observations are found in 
Sect. 3. The hyperfine structure (hfs) components of the nitrogen 
hydrides are discussed in Sect. 4. In Sect. 5 we use three different 
methods to decompose the absorption lines along the sight-lines 
towards both sources in different velocity components, and es- 
timate column densities, N, and relative abundances, X, in each 
component. We also compare the column densities of the nitro- 
gen hydrides, both with each other to investigate possible cor- 
relations, and with other species tracing regions with both low 
and high molecular fractions. Section 6 presents our estimates 
of the ortho-to-para ratio (OPR) of NH3. We end this paper with 
a summary and outlook in Sect. 7. Note that the analysis of the 
background source emissions and absorptions is left for a future 
paper. 

2. Observations and data reduction 

All the reported observations are a part of the more extended 
PRISMAS programme towards G10.6-0.4 and W49N. The ni- 
trogen observations, which took place between March 2 and 
April 18, 2010, are summarised in Table 1. All Herschel identifi- 
cation numbers (OBSID's) are found in Table A. 1 (on-line mate- 
rial). Note that NH Nj = 1 — 1 was observed together with CH 2 
at 946 GHz, and ortho-NH3 J K = 2o - lo was observed together 
with para-NH3 Jk — 2\ - 1 1 (the absorptions from the source are 
97 km s _1 apart). Before the launch of Herschel, no observa- 
tions of these NH 3 transitions, nor the NH + 2 n l/2 J = 3/2 <- 1/2 
transition, had been performed. We do not, however, analyse 
the ortho-NH3 Jk = 2o - lo transition in this paper since this line 
shows absorption only at the background source velocities. 

We used the dual beam switch mode and the wide band 
spectrometer (WBS) with a bandwidth of 4x1 GHz and an ef- 
fective spectral resolution of 1.1 MHz. The corresponding ve- 
locity resolution is about 0.3 kms -1 at the highest frequen- 
cies and 0.6 km s _1 at 572 GHz. In addition, simultaneous ob- 
servations were performed using the high resolution spectrom- 
eter (HRS) with an effective spectral resolution of 0.25 MHz 
(Av~0.1 km s" 1 ) and a bandwidth of 240-340 MHz except for 
the 1214 and 1215 GHz lines which had an effective resolution 
of 0.5 MHz (Av-O.lkms- 1 ) and a bandwidth of 780 MHz in 
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order to cover both lines in the same band. Note that the chan- 
nel separation of the observations is 0.5 MHz in the WBS and 
0.12 MHz in the HRS (0.24 MHz for the 1 215 GHz line). Each 
line was observed with three different frequency settings of the 
local oscillator (LO) corresponding to a change of approximately 
15 km s _1 to determine the sideband origin of the lines. During 
all observations two orthogonal polarisations were used. 

The pointings were centred at a — 19 h 10 m 13 ^2, 

5 = 09° 06' 12" (72000) for W49N, and a = 18 h 10 m 28 s .7, 

6 = -19° 55' 5" (72000) for G10.6-0.4. The recommended 
values for half-power beam width of the telescope are 44'.'2, 
22'.' 1, and 1879 at 480, 960 and 1 120 GHz, respectively. The 
reference beams were located within 3' on either side of the 
source. The total calibration uncertainties are <9% for band 1 
and <13% for band 5, including the sideband gain ratio un- 
certainty which is 3-4% for band lb (ortho-NH 3 ), and 4-6% 
for band 5a (para-NH3). All errors are added in quadrature. 
Detailed information about the HIFI calibration including beam 
efficiency, mixer sideband ratio, pointing, etc., can be found 
on the Herschel internet site 3 . The in-flight performance is 
described by Roelfsema et al. (2012). Since we are interested 
in absorption lines and their relative strength compared to the 
continuum in this paper, the forward efficiency of 96% and 
the main beam-efficiency of 64-76% (between 1 120 and 
480 GHz) do not affect our results and we have therefore used a 
T* K intensity scale throughout this paper. 

The data were processed using the standard Herschel 
Interactive Processing Environment (HIPE, 4 Ott 2010), ver- 
sion 5.1, up to level 2 which provides fully calibrated spectra. 
The data quality is excellent with very low intensity ripples in 
most cases, typically below a few percent of the double side- 
band continuum. For NH + in both sources and both polarisa- 
tions, and 0-NH3 at 572 GHz in the vertical polarisation to- 
wards G10. 6-0.4, we corrected the data from ripples using the 
fitHifiFringe tool in HIPE (in paper I we only used one polarisa- 
tion for these two transitions). The FITS files were then exported 
to the spectral line software package xs 5 which was used in the 
subsequent data reduction. 

The two polarisations are generally in agreement to better 
than 10%. The three LO-tunings are also in very good agree- 
ments without contamination from the image sideband with one 
exception in W49N around the 0-NH3 at 572 GHz line (see 
Sect. 3.1.1). In all other spectra we averaged the three LO- 
tunings and both polarisations. 

For the above line identification work we used the Cologne 
Database for Molecular Spectroscopy 6 (CDMS, Muller et al. 
2005), the Jet Propulsion Laboratory catalogue 7 (JPL, Pickett 
et al. 1998), and Frank Lovas' Spectral Line Atlas of Interstellar 
Molecules 8 (SLAIM). Laboratory measurements of NH transi- 
tions have been performed by e.g. Klaus et al. (1997) and Flores- 
Mijangos et al. (2004), and measurements of o-NH 2 by e.g. 
Muller et al. (1999). A review on ammonia is found in Ho & 
Townes (1983), and an investigation of the hyperfine structure 



3 http://herschel.esac.esa.int/Docs/HIFI/html/ch5.html 

4 HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Centre, and the HIFI, PACS and SPIRE consortia. 

5 Developed by Per Bergman at Onsala Space Observatory, Sweden; 
http : //www . Chalmers . se/rss/oso-en/ observations/data- 
reduction-software 

6 http://www.cdms.de 

7 http://spec.jpl.nasa.gov 

8 http://physics.nist.gov/ 



components of ortho-NH3 at 572 GHz was recently performed 
by Cazzoli et al. (2009). NH + measurements have been made by 
Verhoeve et al. (1986) and Hiibers et al. (2009). 

3. Observational results 

Figures 1 and G. 1 (on-line material) show the averaged double 
sideband WBS spectra of the transitions listed in Table 1 towards 
W49N and G10.6-0.4, respectively (special treatment of 0-NH3 
lo-Oo in W49N is found in Sect. 3.1.1). All detected species 
show very similar absorption patterns over a wide range of ve- 
locities, although W49N has fewer absorbing velocity compo- 
nents in the line-of-sight. No detections of NH + have been found 
in any of the two sources. 

The resulting noise and single sideband (SSB) continuum 
levels, for each transition, are found in Table 1 . Note that since 
HIFI uses double side band (DSB) receivers, the observed con- 
tinuum has to be divided by two to obtain the SSB contin- 
uum. The sideband gain ratio is assumed to be unity through- 
out this paper. This has proven to be a good assumption based 
on PRISMAS observations of saturated absorption lines like HF 
that establish the true zero levels (Neufeld et al. 2010b). 

3.1. Background source emission 

Even though we only analyse the absorption lines in the fore- 
ground clouds in this paper we still have to take the source emis- 
sion into account in the analysis. Several spectra show broad 
emission features from the background sources which extend far 
into the velocity range pertaining to the line-of-sight clouds. 

Ammonia in particular shows strong emission at the back- 
ground source velocities with self-absorption around +13 and 
km s -1 for W49N and G10.6-0.4, respectively. A large part 
of the NH 974 GHz emission at the source velocities in both 
sources is most likely not from NH, but from HCN 7=11-10 
at 974.487 GHz, only 9 MHz above the NH line correspond- 
ing to 2.5 km s _1 . The PRISMAS observations of these sources 
have detected three additional HCN emission lines: 7 = 6-5 
at 532GHz, 7 = 7-6 at 620 GHz, and 7= 10-9 at 886 GHz 
which strengthens our identification. In addition, we detect no 
NH emission in the 946 GHz transition which is a further con- 
firmation that the emission is due to another species. HCN in 
foreground gas is unlikely to produce measurable absorption in 
the 7 = 11-10 transition because the excitation energy of the 
lower state, E\o/k = 234 K, is much too high to be populated at 
the density and temperature of diffuse or translucent gas. 

In the analysis performed in Sect. 5 using three different 
methods, we have removed the broad emission lines that over- 
lap with the absorption in the line-of-sight by means of Gaussian 
fits. 

3.1.1. Emission line contamination 

In W49N, close to the 0-NH3 line at 572 GHz, an emission line 
from the lower sideband, identified as the Jk„,k c = 166,io- I6541 
SO2 line, moves across the ammonia absorption lines in the dif- 
ferent LO-tunings. Thus for the 0-NH3 572 GHz transition in 
W49N we were only able to use the A and C LO-settings in 
Table A. 1, which also had to be cut to remove the parts contam- 
inated by the SO2 line. 

The 0-NH2 absorption towards W49N is contaminated by a 
an emission line in the same sideband from the source. We iden- 
tify the emission line as a blend of three unresolved hfs compo- 
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Fig. 1: Nitrogen hydrides towards W49N: Double sideband WBS 
spectra of NH, ortho-NH^, ortho- and para-NH^, and NH + over 
the LSR velocity range -50 to 85 km s _1 . Quantum numbers are 
found in Table 1 . Note that we only analyse the absorption in the 
velocity range 30-75 km s _1 in this paper and leave the source 
absorption and emission for a future paper. 



nentsofN0 2 n 1/2 J = 9.5 f -» 8.5 f ,F= 10.5-9.5, 9.5-8.5, and 
8.5 - 9.5, at 952.464201 GHz. The NO emission line is removed 
from the 0-NH2 absorption (described in Sect. B, on-line mate- 
rial), and in all following figures in this paper we show 0-NH2 
in W49N with removed NO emission. The NO identification is 
strengthened by the PRISMAS observations of three additional 
NO lines towards W49N shown in Fig. B.2 (on-line material). 
No emission of NO is found towards G10.6-0.4. 




-20 



Velocity [km s 1 ] 

Fig. 2: W31C: o- u NH 3 and o- 15 NH 3 . Both isotopolo gues show 
emission from the source. Absorption in the line-of-sight 
is only seen in the o- 14 NH3 spectra. The emission line at 
vlsr ~41 km s _1 in the o- 15 NH3 spectra, is an SO2 line from the 
upper sideband. 



3.1.2. The 15 NH 3 isotopologue: only in G1 0.6-0.4 

Since our observations also covered the frequency of o- 15 NH3 
/ff = lo-Oo at 572.113 GHz, we have checked the averaged 
spectra for emission and/or absorption features of this isotopo- 
logue. In Figs. 2 and G.2 (on-line material) we show a 5cr emis- 
sion feature in G10.6-0.4 with an amplitude of 40 mK and a 
line-width of about 4 km s _1 , which we identify as o- 15 NH3 

lo-Oo, giving a source velocity of 4 km s _1 . Checking the 

possibility of absorption lines of o- 15 NH3 in the line-of-sight 
towards G10. 6-0.4, we note that there seems to be an absorp- 
tion feature corresponding to the strongest velocity components 
of o- 14 NH3 at ~27-32 km s -1 . This feature is, however, prob- 
ably not a real detection since it would otherwise imply that 
14 N/ 15 N ~ 20. This is about 20 times lower than the 14 N/ 15 N ra- 
tio in the local ISM which is about -450 (Wilson & Rood 1994), 
or 441+5 in the Solar nebula (Marty et al. 201 1). Odin observa- 
tions of absorption towards Sgr B2 show 14 NH 3 / 15 NH 3 > 600 
(Wirstrom et al. 2010). Lower ratios have been found, for in- 
stance, in low mass cores, 334+50 (Lis et al. 2010), and on Earth, 
-270. Adande & Ziurys (201 1) find values of -100-350 using 
CN and HNC in ten dense molecular clouds located at various 
galactic distances from the Galactic Centre, but this is still sev- 
eral times higher than inferred from a possible o- 15 NH3 absorp- 
tion. The 14 N/ 15 N ratio is also difficult to determine and the re- 
sults using cyanides may differ from the ratio inferred from other 
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W49N:NH 974GHz W49N: o-NH 2 953 GHz 





Velocity [km/s] 



Velocity [km/s] 



Fig. 3: NH hfs components: Normalised WBS spectrum to- Fig. 4: Ortho-NH2 hfs components: Normalised WBS spectrum 
wards W49N. Also shown are Gaussian fits with a line width towards W49N. Also shown are Gaussian fits with a line width 
of 1 km s _I of the NH 974 GHz hfs components in the of 1 km s _1 of the 0-NH2 953 GHz hfs components in the 



+39 km s velocity component. 



+39 km s velocity component. 



species such as N2H + , and may not reflect the true 14 N/ 15 N abun- 
dance ratio. More observations are needed to lower the noise in 
order to obtain a real detection of 15 NH3 and to determine the 
true 14 N/ 15 N ratio towards G10.6-0.4. 

In W49N, no o- 15 NH3 detection is made due to blending 
SO2 lines from the upper sideband which are both stronger and 
broader than in G10.6-0.4. In addition, the o- 14 NH3 emission is 
more than two times weaker than in G10. 6-0.4, and the spec- 
trum has a three times higher noise level which does not allow a 
detection of o- I5 NHi emission from the source. 



4. Hyperfine structure components 

Both NH and NH2 have numerous hyperfine structure (hfs) com- 
ponents which require a high resolution spectrometer, such as 
HIFI, to be spectrally resolved. The on-line Tables C.l -C.3 list 
the hfs components of the NH 974 and 946 GHz, and the o-NH 2 
953 GHz transitions which have 21,9, and 30 hfs components, 
respectively. The observed NH 974 GHz and 0-NH2 transitions 
in the +39 km s _1 velocity component towards W49N are shown 
in Figs. 3 and 4 together with models of respective transition, 
including hfs components, using Gaussian opacity profiles and 
a line width of 1 km s _1 . Here, the intensities have been nor- 
malised to the continuum in single sideband as TjJTc-1 assum- 
ing a sideband gain ratio of unity where is the observed in- 
tensity and Tc is the SSB continuum as measured in line-free 
regions in the spectra. The large number of hfs components, in 
addition to many different, overlapping, velocity components in 
the line-of-sight, complicates the analysis substantially. In par- 
ticular, the 0-NH2 hfs components in the line-of-sight absorp- 
tion overlap with the source emission/absorption, thus requiring 
a source model in the analysis. 

The much simpler hyperfine structure of the ortho- and para- 
NH3 ground state rotational transitions has never been resolved 
in space before Herschel. The Odin satellite was able to observe 
the 572 GHz ammonia transition, but the velocity resolution of 
about 0.5 km s _1 was not enough to spectrally resolve the hfs 
components even though the asymmetric line shapes hinted the 
hfs components (e.g. Liseau et al. 2003). Figure 5 shows an ex- 
ample of our observations of the 572 GHz 0-NH3 transition in 
the +39 km s _1 absorption velocity component towards W49N 
with the HRS and a velocity resolution of 0.13 km s _1 . Three 
hfs components are seen at 38.4, 39.4 and 40.0 km s _I although 
the line width of 1.0 km s prevents a detailed check of the 



W49N: 0-NH3 572 GHz high resolution 




Velocity [km/s] 



Fig. 5: Ortho-NH^ hfs components: Normalised high resolution 
spectra (Av = 0. 13 km s -1 ) towards W49N in the +39 km s -1 ve- 
locity component. Three hfs components are seen at 38.4, 39.4 
and 40.0 km s _1 . Also shown are the positions and relative in- 
tensities of the hfs components and their fitted opacity profiles 
with a line width of 1 .0 km s -1 . 



relative intensities. In the Gaussian fit we therefore use 0.2, 1 
and 0.6 as relative intensities, with corresponding relative off- 
sets at -1.05, 0.0, and 0.52 km s" 1 found by Cazzoli et al. (2009) 
(Table C.4, on-line material). The positions and relative inten- 
sities of the hfs components and the fitted opacity profiles are 
also shown in Fig. 5. The six hfs components in the para-NH? 
line at 1215 GHz are even more closely spaced and not resolved 
(Table C.5). 

The WBS ammonia observations do not spectrally resolve 
the NH3 hfs components (Av = 0.6 km s _1 at 572 GHz and 
Av = 0.3 km s _1 at 1215 GHz), but we still take the hfs compo- 
nents into account in the following NH3 modelling since they 
produce a slightly asymmetric line shape and a systematic line 
broadening. 



5. Analysis: Velocity decomposition, column 
densities, and abundances 

In paper I we estimated the abundances towards G10.6-0.4 sim- 
ply by evaluating the integrated opacity of each line over the ve- 
locity range ll-54kms~ 1 . This approach gave a first estimate 
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of the average abundance in the full line-of-sight, but did not 
take into account possible abundance variations in the different 
velocity components. In addition, no source model was included, 
and thus the blend of the NH and 0-NH2 hfs components from 
the foreground absorptions with the source absorption was ne- 
glected. 

In this paper we use three different approaches to decompose 
the absorption into separate velocity components, and to ob- 
tain column densities and abundances in each component. Each 
method has its own strengths and weaknesses, and differences 
in the results can be considered an estimate of the errors of the 
methods. If all three methods point to the same result, it will be 
considered as robust. We note, however, that the uncertainties of 
the data are higher towards W49N than for G10.6-0.4, since we 
have removed a SO2 line from the upper sideband in the ortho- 
NH 3 absorption, and an NO line in the o-NH 2 absorption. 

In this work we also only estimate N and X for the ortho- 
symmetries of NH2 and NH3 since our own measurements of 
the ammonia OPR in Sect. 6 do not point to the high tempera- 
ture limit, but instead suggest an OPR lower than unity, which 
cannot be straightforwardly explained. If this unexpected low 
ammonia OPR is true, then the processes that affect the ammo- 
nia OPR perhaps also could affect the OPR of NH2 in diffuse 
clouds about which we have no information. And since we want 
to compare relative abundances, we choose to focus on the ortho- 
symmetries since both are observed in NH2 and NH3, and also 
since the 0-NH3 observations have higher S/N than the P-NH3. 

To obtain integrated opacities of each transition, in each 
velocity component, the first method uses Gaussian fitting 
(Sect. 5.2), while the second one uses the observed spectra of 
ortho-NH3 and CH as templates for other species (Sect. 5.3). 
Both methods calculate the opacity as r=-ln(7A/7c - 1). The 
derived integrated opacities are then used to estimate column 
densities by means of the non-LTE (Local Thermodynamic 
Equilibrium) RADEX code (Sect. 5. 1). As a third method, we have 
used XCLASS (Sect. 5.4) which, in contrast to RADEX, assumes 
fixed excitation temperatures, to obtain column densities in each 
velocity component directly. 

In order to be able to use our results as comparisons with 
chemical models, we need to derive the abundance of the species 
with respect to the total amount of gas. The difficulty is to obtain 
reliable estimates of N(H) and/or MFb). In paper I, we com- 
pared the nitrogen hydrides column densities to the total amount 
of hydrogen in the line-of-sight, N H = N(H)+2jV(H 2 ). This can 
not be done in the present work since the neutral hydrogen ab- 
sorption shows too broad and overlapping absorption compo- 
nents over much larger velocity ranges as compared to the nitro- 
gen hydrides. In order to estimate the amount of molecular hy- 
drogen in each (narrow) velocity component, we therefore con- 
sider the comparison with CH (Gerin et al. 2010b) as an abun- 
dance determination. This can be made using the CH/H2 corre- 
lation (CH/H 2 = 3.5xl0- 8 ) observed by Sheffer et al. (2008) in 
the solar neighbourhood diffuse medium. This correlation is as- 
sumed to be valid in regions which are dominated by ultraviolet 
radiation where N(CH) < 2xl0 14 cm -2 . In this way we can use 
our own PRISMAS observations of CH and measure abundance 
ratios obtained with the same instrument. 

The physical properties of the bulk of the absorbing 
gas in our sight-lines are (Gerin et al. 2010b), according 
to the definitions of Snow & McCall (2006), a mixture of 
diffuse (Nu < 500 cm 4 and A V <1) and translucent gas 
(500 < N H <5 000 cm" 3 and a shielding 1 < A v <5). And since 
the comparison later in this section of the nitrogen hydrides with 
other species tracing both high and low density gas, points to an 



existence in the denser parts of the gas, molecular hydrogen is 
most likely the dominant form of hydrogen in these components. 

Assuming that the nitrogen hydrides co-exist with HF and 
H2O in the line-of-sight material we use the fact that the satu- 
rated absorptions of the latter two species reach the zero level, 
for a sideband ratio of unity, to support our assumption that 
the foreground absorbing material completely covers the contin- 
uum within the beam. Comparison spectra of 0-NH3 and o-H 2 
are shown in Figs. 8-9, and with HF, ortho- and para-H20 in 
Figs. G.9-G.10 (on-line material). 



Table 2: Conversion factors'" . 



Species x 


Transition 


T 

1 ex 


N(x)(J rdv= 1.0km s" 1 ) 






(GHz) 


(K) 


(cm- 2 ) 




NH 


946.476 


3.9 


3.65x10" 


a 




974.478 


4.0 


7.47xl0 12 


a 


0-NH2 


952.578 


3.9 


3.62xl0 12 


b 


0-NH3 


572.498 


2.8 


3.92xl0 12 


b 


p-NH 3 


1 215.245 


4.9 


1.65xl0 13 


b 


NH + 


1 012.524 


4.2 


6.48xl0 12 




CH 


532.724 


2.9 


2.28xl0 13 


b 


HNC 


90.664 


2.8 


1.80xl0 12 


b 


CN 


113.169 


2.8 


1.90xl0 13 


b 



Notes. <a) The excitation temperatures and conversion factors from in- 
tegrated opacity, Jrdv = 1.0 km s -1 , to column densities, N(x), are cal- 
culated using RADEX with T k = 30 K and «(H 2 ) = 500 cnT 3 which cor- 
responds to n H <; 10 3 cm 3 . (A) Integrated over all hfs components. 



5.1. RADEX 

Observations of absorption lines can provide accurate determi- 
nations of molecular column densities if the observed transi- 
tions trace all of the populated states. In the case of NH 3 we 
observe two transitions of the ortho (A-symmetry) form, aris- 
ing in the Jk = % and lo states, and one transition of the para 
(E-symmetry) form arising in the lower li rotation-inversion 
level. Although the upper 1 1 inversion level is less metastable 
than the lower one, it may still be significantly populated. The 
metastable 22 para levels and 33 ortho levels may also have 
significant populations in diffuse molecular clouds at low den- 
sity, n H < 10 4 cm 3 , owing to infrared pumping and formation. 
Because we are directly sensitive to three lower states and lack 
direct information about the excitation temperatures of the ob- 
served transitions, we use the non-equilibrium RADEX code (van 
der Tak et al. 2007) to relate the total column density of ammonia 
to the observed integrated optical depths. The models provide in- 
tegrated opacities of both observed and unobserved excited lev- 
els and quantify possible effects of stimulated emission, chem- 
istry, and electron collisions. Results for one reference model are 
summarised in Table 2, where conversion factors for all nitrogen 
hydrides and comparison species are listed in terms of the total 
column density N(x) of molecule x that is needed to achieve an 
integrated optical depth J rdv -1.0 km s 1 in the specified tran- 
sition. In this model we treat ortho- and para-symmetries sepa- 
rately. The integrated optical depth is a sum over all hyperfine 
structure in the transition. 

We have assumed diffuse molecular cloud conditions with a 
kinetic temperature of 30 K and a density of molecular hydrogen 
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n(H2) = 500 cm" 3 , which corresponds to n# <: 10 3 cm" 3 . The re- 
sulting conversion factors for the rotational ground state transi- 
tions are not very sensitive to density, temperature or electron 
collisions. The critical densities of the fundamental rotational 
transitions of NH, NH2, and NH3 are high, n cr j t ~10 8 cm" 3 . This 
means the upper energy levels of these transitions must be ex- 
cited almost entirely radiatively at the low densities in diffuse 
gas. We find, however, that electron collisions can be responsi- 
ble for excitation temperatures of ~4 K in the cm-wave inversion 
transitions rather than < 3 K that would be found in conventional 
excitation analyses at low densities, n(H2) < 10 3 cm" 3 . Where 
collision rates are unknown, we have made guesses scaled in 
proportion to radiative line strengths. The background radiation 
field includes both the 2.725 K cosmic microwave background 
and a model of the Galactic infrared radiation in the solar neigh- 
borhood. The resulting excitation temperatures of the observed 
sub-millimetre transitions are typically 3-4 K, which are small 
enough compared to hv/k that no correction for emission is re- 
quired. Note that the lower state of the 572 GHz transition con- 
tains 95% of the ortho molecules while the lower state of the 
1215.2 GHz transition occupies only 43% of the para molecule. 

5.2. Method I: Multi-Gaussian fit of the nitrogen hydrides 
simultaneously 

We have modelled the observed spectra of NH, 0-NH2, and 
0-NH3 using Gaussian optical depth profiles. These were gener- 
ated for each hfs component of each species in separate velocity 
intervals, and were made to fit the observations under the condi- 
tion that the LSR velocity and width in each velocity component 
must be the same for all molecules. This can be made assuming 
that all species co-exist in the same velocity ranges and there- 
fore show absorption at the same velocities. This assumption is 
supported by the striking similarities of the strongest NH, 0-NH2 
and 0-NH3 absorptions, despite the complicated hyperfine struc- 
ture patterns of NH and 0-NH2, shown in normalised compari- 
son spectra in Fig. G.9 (on-line material). We note that a good fit 
therefore is a good indicator if this assumption is valid. 

The observed line profiles are thus modelled by Gaussian 
components of both the hfs components and all velocity compo- 
nents according to 



= ex P 



~ 2 X T0(0 6XP[ ~ 4 ln(2) 

i=l 7=1 



V - V (0 - Vms(7) 



Av 



where N v is the number of the modelled velocity components 
(same for all species), Nut* is the number of hfs components 
(which differs between species), Thf s is the theoretical relative 
line strength of each hfs component, To is the opacity of the 
strongest hfs component in each velocity component, vo is the 
LSR velocity of each velocity component, Vhf s is the relative ve- 
locity offset of the hfs components with respect to the strongest, 
and Av is the FWHM (full width of half maximum). We assume 
that Thfs = A u i x g„/(A„;(main) x g„(main)) for the hfs compo- 
nents of NH and 0-NH2, since the many overlapping velocity 
components prevent a check of the relative intensities of these 
components. 

We have used both ortho- and para-NH3, which are not sig- 
nificantly complicated by hfs splitting, to determine the mini- 
mum number of velocity components needed, and also to set-up 
reasonable initial guesses for the line properties necessary for 
the fitting procedure to converge at all. The fitting was done for 
each velocity component at a time in a loop until all fits had con- 
verged. The fits are most likely not unique. Provided that the fits 



are good and that we use the same velocity parameters for all 
species, this is, however, not considered important as we do not 
ascribe any physical meaning to the Gaussian components. 

The results from Method I are peak opacities, line widths 
and centre velocities for each velocity component and species. 
The integrated opacity in each velocity component was assumed 
to be well represented by ^Tdv = 1.06AvT pea k- The non- 
detections of NH + are used to put upper limits on its integrated 
opacities (3cr). We then used the conversion factors in Table 2 to 
calculate the column densities tabulated in Table 3. Typical er- 
rors in the G10.6-0.4 resulting column densities for NH, 0-NH2 
and 0-NH3 are between 7 and 15%. In the vlsr = 45 km s" 1 com- 
ponent, the errors are 22% for all three species, and 23% in the 
vlsr = 39 km s" 1 component for NH. The errors in W49N varies 
between 8 to 1 1% for all three species, except in the 33 km s" 1 
component where the uncertainties are 37, 26 and 31% for NH, 
0-NH2 and 0-NH3, respectively. The errors include calibration 
uncertainties (see Sect. 6) as well as uncertainties from the fit- 
ting procedure. 

In addition to the nitrogen hydrides, we also modelled the 
PRISMAS CH 532 GHz transition, and the CN 113 GHz and 
HNC 91 GHz transitions observed with the IRAM 30 m antenna 
(Godard et al. 2010). These species were not included in the fit- 
ting procedure described above. Instead we have used the result- 
ing velocity parameters from the fits of the nitrogen hydrides 
(vlsr an d li ne widths) since we want to investigate possible cor- 
relations and make abundance determinations only in the same 
parts of velocity space as the nitrogen hydrides. The resulting 
CH column densities are found in Table 3, and the CN and HNC 
column densities in the on-line Table F. 1 . 

Figures 6 and 7 show the fits and residuals of the nitrogen 
hydrides together with CH. The resulting fits are well reproduc- 
ing the observations for all species except for CH suggesting that 
CH exist in a more widely spread gas than the nitrogen hydrides. 

The CH column densities are then used together with the 
relation [CH]/[H 2 ] = 3.5x1 0" 8 (Sheffer et al. 2008) to estimate 
the abundances listed in Table 3. 



5.3. Method II: ortho-NHj, and CH as templates for other 
species 

This method directly compares the opacity line profile of the 
rtoh'- saturated 0-NH3 line at 572 GHz with the line profiles 
of NH and 0-NH2, convolved with respective hyperfine struc- 
ture. This assumes that all species co-exist in the same velocity 
ranges, and that the opacity ratios depend only on their column 
density ratios, which may vary from one velocity component to 
another. No opacity profile is assumed for a given absorbing 
cloud so, unlike what is done by methods I and III, we do not 
attempt to decompose the velocity structure of the template line 
into a number of Gaussian curves. We perform, instead, simple 
cuts of the template profile into large velocity bins. 

The hfs of the 0-NH3 line is not resolved in the WBS data 
and introduces only a small broadening of the lines that is dealt 
with by smoothing the NH and 0-NH2 spectra. Before compar- 
ing the opacity profiles, we mask out those portions of the spec- 
tra where a simple scaling between template and target is not ex- 
pected, such as the velocity range associated to the background 
sources, and remove the emission line wing by Gaussian fitting. 
We note that the removal of the line wing introduces very large 
uncertainties towards W49N around 33 km s" 1 where a weak ab- 
sorption component is visible. This method therefore disregards 
the results in this velocity component. 
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Table 3: Method I results' 1 . 



W49N 

vlsr Av AT(NH) JV(oNH 2 ) JV(oNH 3 ) A'(NH + ) Af(CH) JL §| *(NH)* X(oNH 2 )* X(oNH 3 ) fc X(NH + ) 6 



(kms- 1 ) 


(kms-') 


(cm 2 ) 


(cm 2 ) 


(cm 2 ) 


(cm" 2 ) 


(cm" 2 ) 














33.3 


2.1 


1.5el2 


2.1el2 


4.6ell 


<3.3ell 


1.6el3 


3.3 


4.6 


3.3e-9 


4.6e-9 


1.0e-9 


<7.2e-10 


39.5 


1.1 


1.0el3 


7.3el2 


4.3el2 


<2.4ell 


2.3el3 


2.3 


1.7 


1.5e-8 


l.le-8 


6.5e-9 


<3.7e-10 


59.2 


2.9 


9.2el2 


4.1el2 


2.1el2 


<7.5ell 


6.0el3 


4.4 


2.0 


5.4e-9 


2.4e-9 


1.2e-9 


<4.4e-10 


62.7 


2.3 


8.0el2 


3.9el2 


2.2el2 


<6.3ell 


5.1el3 


3.6 


1.8 


5.5e-9 


2.7e-9 


1.5e-9 


<4.3e-10 


Total: 




2.9el3 


1.7el3 


9.1el2 


<2.0el2 


1.5el4 


3.2 


1.9 


6.7e-9 


4.1e-9 


2.1e-9 


<4.6e-10 


Mean: 




7.2el2 


4.4el2 


2.3el2 


<4.9ell 


3.8el3 


3.2 


1.9 


6.7e-9 


4.1e-9 


2.1e-9 


<4.6e-10 


Median: 




8.6el2 


4.0el2 


2.2el2 


<4.8ell 


3.7el3 


4.0 


1.9 


8.1e-9 


3.8e-9 


2.0e-9 


<4.5e-10 



G10.6-0.4 

v LS R Av AT(NH) JV(oNH 2 ) JV(oNH 3 ) N(NH + ) N(CH) Jg. ||l X (NH) b X(oNH 2 ) b X(oNHi) b X(NH + ) b 

(kms -1 ) (kms -1 ) (cm -2 ) (cirr 2 ) (cm -2 ) (cm -2 ) (cm -2 ) 



16.2 


1.7 


2.0el3 


9.6el2 


4.2el2 


<4.7ell 


5.8el3 


4.8 


2.3 


1.2e-8 


5.8e-9 


2.5e-9 


<2.8e-10 


18.8 


1.5 


6.1el2 


3.1el2 


8.7ell 


<4.1ell 


2.6el3 


7.0 


3.6 


8.2e-9 


4.2e-9 


1.2e-9 


<5.5e-10 


22.1 


4.4 


1.6el3 


6.9el2 


1.8el2 


<1.2el2 


7.7el3 


8.9 


3.8 


7.3e-9 


3.1e-9 


8.2e-10 


<5.5e-10 


22.8 


1.0 


1.0el3 


3.6el2 


1.5el2 


<2.7ell 


7.1el2 


6.7 


2.4 


4.9e-8 


1.8e-8 


7.4e-9 


<1.3e-9 


24.8 


3.1 


6.2el2 


4.3el2 


1.4el2 


<8.4ell 


3.7el3 


4.4 


3.1 


5.9e-9 


4.1e-9 


1.3e-9 


<8.0e-10 


27.8 


1.9 


2.7el3 


1.2el3 


7.3el2 


<5.0ell 


7.4el3 


3.7 


1.6 


1.3e-8 


5.7e-9 


3.5e-9 


<2.4e-10 


29.9 


1.7 


2.9el3 


1.3el3 


7.9el2 


<4.6ell 


5.6el3 


3.7 


1.6 


1.8e-8 


8.1e-9 


4.9e-9 


<2.9e-10 


32.1 


1.5 


9.0el2 


3.5el2 


1.3el2 


<4.1ell 


3.5el3 


6.9 


2.7 


9.0e-9 


3.5e-9 


1.3e-9 


<4.1e-10 


36.1 


4.0 


2.0el3 


6.5el2 


3.2el2 


<l.lel2 


7.8el3 


6.3 


2.0 


9.0e-9 


2.9e-9 


1.4e-9 


<4.9e-10 


39.0 


1.2 


l.lel3 


3.0el2 


4.7ell 


<3.1ell 


3.3el3 


23 


6.4 


1.2e-8 


3.2e-9 


5.0e-10 


<3.3e-10 


40.9 


1.7 


3.0el3 


l.lel3 


4.1el2 


<4.7ell 


5.4el3 


7.3 


2.7 


1.9e-8 


7.1e-9 


2.7e-9 


<3.0e-10 


45.0 


1.4 


3.0el2 


l.lel2 


5.9ell 


<3.6ell 


l.lel3 


5.1 


1.9 


9.6e-9 


3.5e-9 


1.9e-9 


<l.le-9 


Total: 




1.9el4 


7.8el3 


3.5el3 


<6.8el2 


5.5el4 


5.4 


2.2 


1.2e-8 


5.0e-9 


2.2e-9 


<4.4e-10 


Mean: 




1.6el3 


6.5el2 


2.9el2 


<5.7ell 


4.6el3 


5.4 


2.2 


1.2e-8 


5.0e-9 


2.2e-9 


<4.4e-10 


Median: 




1.4el3 


5.4el2 


1.7el2 


<4.7ell 


4.6el3 


8.2 


3.3 


1.0e-8 


4.2e-9 


1.3e-9 


<3.6e-10 



Notes. (a) LSR velocity, vlsr, line width, Av, column densities, N(x), and abundances, X(x) = 7V(x)/./V(H 2 ), of each velocity component, including 
3a- upper limits of N(NH + ). Note that we estimate N and X only for the ortho-symmetry forms of NH 2 and NH 3 . {b) Using [CH]/[H 2 ] = 3.5xl0" 8 
(Sheffer et al. 2008). 



The velocity bins into which the 0-NH3 opacity profile is 
split are chosen to be significantly larger than the original resolu- 
tion, but narrow enough to separate the most prominent velocity 
features (typically a few km s _1 , see Table D.l). Each velocity 
bin of the 0-NH3 template is then convolved, channel by chan- 
nel, with the NH and 0-NH2 hfs components. We thus obtain an 
intermediate opacity model of how the NH and 0-NH2 spectra 
would appear if absorption were limited to the selected velocity 
bin. The observed 0-NH3 opacity is then fitted by the IDL Least 
Squares Fitting routine MPFIT (Markwardt 2009) with a linear 
combination of the intermediate models. The output of the fits 
consists of the opacity ratios, between the 0-NH3 and the mod- 
elled NH and 0-NH2 spectra, in the various velocity bins. 

To obtain column densities in each velocity bin, we sum the 
opacities of all channels of the template transition, and then ap- 
ply the RADEX conversion factors from Table 2. The resulting 
abundance ratios are finally used to obtain column densities of 
the other species. Typical errors (including both calibration and 
fitting uncertainties) are 5 - 15% except for the 43.5-50 km s _1 
velocity bin in G10.6-0.4 and the 50-57 km sT l velocity bin in 
W49N which both have an uncertainty of approximately 25%. 

The main advantage of this method is that a small number 
of free parameters is needed to obtain a reasonable fit since it 
uses, as much as possible, the information carried by the tem- 
plate spectrum. The model also avoids the problem of non- 



uniqueness of Gaussian decompositions of the absorption pro- 
files and allows a straightforward comparison of column den- 
sities evaluated for exactly the same velocity intervals. On the 
other hand, such rigidity in the definition of the velocity bins 
does not allow us to reliably retrieve the column density ratios 
if the absorption comes from very different gas volumes result- 
ing in significant velocity shifts of the absorption features be- 
tween the two species. We note that the fit quality, estimated with 
J |residuals| x dv / Jt x dv, is a good indicator of such cases. 
Note also that residual baseline structures, like those produced 
by standing waves or by an imperfect removal of ammonia emis- 
sion line wings, may introduce artificial discrepancies between 
the template spectra and the modelled species. 

In addition to the above modelling, we use the 0-NH3 as a 
template for CN and HNC (Godard et al. 2010), and also use the 
deconvolved CH spectrum at 532.724 GHz (Gerin et al. 2010b) 
as a template to model all three nitrogen hydrides, in order to 
obtain an estimate of the abundance with respect to molecular 
hydrogen. 

The resulting column densities and relative abundances in 
different velocity bins are given in Table D.l, except for the CN 
and HNC results which are found in on-line Table F.2. 

The model fits compared to the normalised spectra in both 
sources are shown in Figs. G.3-G.6 (on-line material). The fit 
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Fig. 6: W49./V: Fits of Method I and the SSB normalised spectra 
of the nitrogen hydrides and CH. The black lines are the observa- 
tions, and the coloured lines are the model fits. In the bottom the 
residuals are plotted on top of each other with respective colour. 
The three transitions used in Method I to determine the vlsr and 
line widths of the velocity components are marked in bold. 



qualities of the comparisons between the 0-NH3, NH and 0-NH2 
are excellent, except in the narrow and deep +39 km s com- 
ponent towards W49N which was very difficult to model. The 
modelling of the nitrogen hydrides using CH as a template show 
on the other hand rather low fit qualities, suggesting that the as- 
sumption that the nitrogen hydrides and CH only come from the 
same gas is less justified, also supported by Method I. 



Fig. 7: G10.6-0.4: Fits of Method I and the SSB normalised 
spectra of the nitrogen hydrides and CH. Same notation as in 
Fig. 6. 



5.4. Method III: XCLASS 

Finally, we have generated synthetic spectra separately for NH, 
0-NH2, 0-NH3, CH, CN, and HNC using the software XCLASS 9 
(created by P. Schilke, and described by Comito et al. 2005, 
and references therein) by considering for each species all hfs 
components and velocity components and the assumption of 
a fixed excitation temperature. The synthetic spectra generated 
with XCLASS were fitted to the observed spectra using MAGIX, 
an iterating engine that allows automatic minimization to con- 
strain model parameters. XCLASS accesses the CDMS and JPL 
molecular databases, and models each molecule with the follow - 



9 https : //www. astro . uni-koeln.de/proj ects/schilke/ 
XCLASS 
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ing free parameters: source size, temperature, column density, 
line width and velocity offset relative to the systemic velocity of 
the source, and derives the column density corresponding to the 
different velocity components detected in absorption and emis- 
sion in the observed spectra. The source size refers to the relative 
size of the absorbing cloud vs. the continuum source, which is 
assumed to be fully covered within the beam. We have assumed a 
Gaussian profile for each hfs component, and an excitation tem- 
perature of 4 K. Since Method III explicitly assumes a fixed ex- 
citation temperature for the estimate of column densities, it is 
significantly different from methods I and II. 

Figures G.7 and G.8, in the on-line material, show the 
model fits together with the observed SSB normalised spectra. 
Resulting column densities and abundance ratios of the nitrogen 
hydrides and CH are found in the on-line Table E. 1, and the CN 
and HNC column densities in on-line Table F. 1 . 



5.5. Comparison of results 

The results from all three methods agree very well, especially 
considering that the velocity bins used are not the same. 

The total N(CH) from all methods are lower than found 
by Gerin et al. (2010b), which is expected since we only use 
parts of the CH absorptions. Comparisons of CH in different 
velocity bins give much more similar results in both sources. 
Since we use CH as a tracer of H2, we also compare our esti- 
mate of N(U 2 ) with Nyman & Millar (1989) towards W49N: for 
the +39 km s _1 component they suggest Af(H 2 ) = 3xl0 21 cm -2 
and they also quote jV( 13 CO) = 7-30xl0 15 cirT 2 which gives 
N(H 2 ) =7-30xl0 21 cirT 2 using a standard conversion fac- 
tor of 10~ 6 . They also suggest column densities towards 
W49N between SO-TOkms" 1 : AT(H 2 ) = 5xl0 21 crrT 2 and 
JV(H 2 )=5.5-23xl0 21 cm" 2 from 13 CO observations. All these 
suggestions for N(H 2 ) are several times higher than inferred 
from our CH observations. This is, however, not surprising be- 
cause the 13 CO conversion factor is rather uncertain at such 
low column densities. Plume et al. (2004) also obtain higher 
JV(H 2 ) than we do towards W49N. The total amount of iV(H 2 ) 
in the sight-lines towards our sources, can also be compared 
with estimates from the K-band extinction (Marshall et al. 
2006): JV(H 2 ) = 1.3xl0 22 and7.5xl0 21 crrr 2 towards W49N and 
G10. 6-0.4, respectively. 

The resulting NH mean abundances of all three methods in 
all velocity components are 5.5 xlO" 9 and l.lxlO" 8 in W49N 
and G10. 6-0.4, respectively. This can be compared to the av- 
erage Af(NH)/jV(H 2 ) value of 3xl0~ 9 in diffuse and translucent 
sight-lines found by Weselak et al. (2009). The mean values of 
o-NH 2 are 3.1xl0~ 9 and 4.5xl0~ 9 in W49N and G10.6-0.4, re- 
spectively, and 1.5xl0~ 9 and 1.9xl0~ 9 foro-NH 3 . 

Our upper limits on the NH + abundances, relative to molec- 
ular hydrogen, are similar in both sources with mean value 
Af(NH + )/JV(H 2 ) <4xl0~ 1() . This is orders of magnitudes lower 
than previous findings from ultraviolet observations towards 
pOph (Snow 1979; de Almeida & Singh 1982), but still much 
higher than the predictions from the chemical models presented 
in Paper I: approximately 10 -13 - 10~ 14 . The upper limits of the 
NH + abundance compared to NH are <2-14% in the differ- 
ent velocity components (except in the uncertain velocity fea- 
ture around +33 km s _1 in W49N), with a mean ratio of <6%. 
This is in contrast to the behaviour of CH + and marginally with 
that of OH + (Wyrowski et al. 2010; Gerin et al. 2010a) with re- 
spect to the corresponding neutrals. The CH + radical reaches 
comparable column densities to those of CH (e.g. Crane et al. 



1995), and JV(OH + ) is a factor 30 below N(OH) in the visible 
data (Krelowski et al. 2010). 

Since our different approaches to deconvolve the velocity 
components for the observed species agree well, we estimate 
that the uncertainties in our derived absolute column densities 
are < 20-50%. The column densities for CH have larger un- 
certainties, but we estimate that the absolute N(CH) results are 
correct within a factor of approximately two. Furthermore, the 
scatter in the CH to H 2 relationship is estimated by Sheffer et al. 
(2008) to a factor of 1.6. In summary, we believe that the abun- 
dance determinations from CH are correct within a factor of a 
few. The abundance ratios of the nitrogen hydrides, relative to 
each other, are on the other hand more accurately determined. 

The results of all three methods confirm our conclusion of 
the analysis of G10. 6-0.4 in paper I that NH is more abundant 
than NH 2 and NH 3 by a factor of approximately 1 .5 - 3.5, assum- 
ing the high temperature ortho-to-para limits of three and one, 
respectively. Note that the mean abundances of 0-NH3 are sim- 
ilar towards both sources, ~2xl0~ 9 , in contrast to o-NH 2 , and, 
in particular NH, which have higher mean abundances towards 
G10.6-0.4 than towards W49N. 

5.6. Comparison and correlations of species 

Figures 8 and 9 show 0-NH3 and five SSB normalised com- 
parison spectra, and in the on-line material Figs. G.9- G.10 
we show additional comparisons of the nitrogen hydrides with 
the deconvolved CH absorption, the HI 21 cm line, H 2 + , OH + , 
HNC, CN, ortho- and para-H 2 0, HF, and HCO + . 

The CH spectra show absorption over a much wider range of 
velocities in both sources than the nitrogen hydrides. Ammonia 
largely follows the CH absorptions, but there are also large 
differences in some parts of the spectra where there is no or 
very little absorption of NH3, while CH shows a much broader 
and stronger absorption, for instance at 35-38, 41-55 and 
65 -70 km s~' towards W49N. Ammonia seems to trace CH 
slightly better towards G10.6-0.4 than W49N. This comparison 
suggests that CH exists in both relatively low and high density 
gas, while the nitrogen hydrides only exist in the parts of the 
interstellar gas with a relatively high density. 

These differences are even more pronounced when compar- 
ing with neutral hydrogen as observed by the VLA A21 cm ab- 
sorption by Fish et al. (2003), which shows absorption over a 
very wide range of velocities with a resolution of 2.5 km s _1 . The 
much more extended velocity space coverage of the HI absorp- 
tion is expected since not all foreground clouds have molecular 
gas. 

The ammonia absorption also follow similar trends as HCO + 
J = l<-0 (Godard et al. 2010), and H 2 0, which is known to 
trace clouds of high molecular fraction (Neufeld et al. 2010b; 
Sonnentruckeret al. 2010). HCO + and H 2 seem to have a rather 
constant abundance ratio in these sight-lines, in contrast to their 
abundances with respect to ammonia. The comparison of am- 
monia to H 2 + and OH + , which mostly reside in lower density 
gas containing considerable amounts of atomic hydrogen (Gerin 
et al. 2010a; Neufeld et al. 2010a), shows no similarities. 

When we compare ammonia with CN and HNC we 
find very similar absorption patterns. The CN 7=1<— 
F = 1/2-3/2, 1/2-1/2 single hyperfine component and the HNC 
J = 1 <— line were observed with the IRAM 30 m antenna 
(Godard et al. 2010). Note, that HNC has three hfs components 
that, however, lie very close in velocities (0.2 and 0.5 km s~'), 
which mostly leads to a broadening of the absorption features. 
These species mainly reside in denser gas than CH, and both 
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Fig. 8: Line-of-sight absorptions towards W49N: Normalised 
(SSB) spectra of ortho-ammonia and species tracing both high 
and low molecular fractions. Strong similarities are seen be- 
tween ammonia, CN, HNC and water, but not with H20 + and 
HI tracing the atomic diffuse gas. 



are closely connected to the NH and NH2 chemistry (e.g. Nejad 
et al. 1990; Liszt & Lucas 2001). 

5.6.1. Column density correlations 

In order to quantitatively examine abundance correlations of the 
nitrogen hydrides, CH, HNC and CN, we show column den- 
sity plots in Figs. 10-11. We have here plotted the results from 
Method I. The parameters of the linear least square fits to the 
data are found in the figures in addition to the correlation coef- 
ficient R. In the on-line material we also show column density 
plots with all three methods (Figs. H.l -H.4). 

Column density plots of the nitrogen hydrides are shown 
on the left hand side in Fig. 10. We note that the scatter is 
larger for jV(NH) vs. _/V(o-NH 3 ) than compared to _/V(o-NH 2 ) vs. 
iV(o-NH3) which show a rather tight correlation. This may indi- 
cate that NH does not entirely exist in the same gas as 0-NH2 




30 40 
v LSR [km s- 1 ] 



Fig. 9: Line-of-sight absorptions towards G10.6—0.4: Notation 
as in Fig. 8. 



and 0-NH3. There is also a possibility that the NH and 0-NH3 
correlation may only be valid for Af(o-NH3) < 5xl0 12 cm" 2 . 
The A^(NH) towards G10.6-0.4 appear to increase up to a max- 
imum of ~3xl0 13 cm" 2 . Also Mo-M-L;) show a slight ten- 
dency of this behaviour towards G10.6-0.4 with a maximum 
of ~1.2xl0 13 cm" 2 . This possible "chemical saturation" corre- 
sponds to a few Ay (estimated from our CH measurements and 
A^h ~ 1.7 x 10 21 cm" 2 mag" 1 x Ay) . This may be explained by 
a more efficient ammonia production at higher column densities, 
or that NH, and 0-NH2, are somehow consumed in the ammonia 
formation. 

On the right hand side in Fig. 10 we show column density 
plots of the nitrogen hydrides vs. CH. The spread of the data 
points reflects the difficulty of finding a velocity structure com- 
mon to hydrides and CH. 

In Fig. 1 1 we show column density plots of the nitrogen hy- 
drides vs. CN on the left and vs. HNC on the right. All three 
nitrogen hydrides show linear correlations with both CN and 
HNC. In dark cloud chemistry, CN and HNC are closely related 
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Fig. 10: Column density comparison plots. On the left hand side: NH, 0-NH2 and 0-NH3. On the right hand side: nitrogen hydrides 
vs. CH. The blue lines show linear least square fits to all data points towards both sources. The correlation coefficients and the linear 
fit equations are given in respective figure. 



to NH and NH 2 through the reactions of C + NH 2 -> HNC + H, 
N + CH -> CN + H, and C + NH -> CN + H (Nejad et al. 1990). 
Also Weselak et al. (2009) found a slightly better correlation be- 
tween NH and CN than compared to CH, and no correlations 
with species such as CH + . 



6. Ortho-to-para ratio of NH 3 

There are two ways to produce molecules in interstellar space: 
in the gas-phase or on grain surfaces. If ammonia is formed in 
the gas-phase, by highly exoergic processes, the ammonia ortho- 
to-para ratio has been expected to be very close to the statistical 
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Fig. 11: Column density comparison plots. On the left hand side we show comparisons of the nitrogen hydrides with CN, and on the 
right hand side comparisons with HNC. Notation as in Fig. 10. 



equilibrium value of 1.0 (the spin statistical value is 4-2 for 
ortho-para, but the number of para states is on the other hand 
almost a factor of two larger). If ammonia is formed at temper- 
atures lower than 40 K on cold dust grains, and then desorbed 
when the grains are heated above 100 K, the OPR may differ 
from unity since the lowest ortho level is 22 K below the lowest 
para level. If no conversion processes between the two symme- 



tries exist, the OPR of ammonia is expected to increase above 
unity at low formation temperatures. 

Prior to Herschel, ortho-to-para ratios have been derived 
from measurements of inversion transitions, the (1,1) and (2, 2) 
transitions involving para states, and the (3, 3) and (6, 6) tran- 
sitions probing highly excited, metastable ortho states. For ex- 
ample, Umemoto et al. (1999) have found OPR= 1.3-1.7 in 
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Fig. 12: W49N. (Upper) WBS spectra of l -0 ortho-NH 3 and 
2 1 — 1 1 para-NH 3 normalised to SSB continuum. (Lower) The 
optical depth ratios are shown for absorptions larger than 3<x as 
a function of LSR velocity. The column density ratios (OPR), 
estimated with RADEX, are given by the right hand y-axis. The 
horisontal dashed line marks an ortho-to-para optical depth ratio 
of 4.2 corresponding to a column density ratio of unity. 



the L 1157 outflow from the observation of six inversion lines 
(J,K) = (1, 1) to (6,6). Similarly, Nagayama et al. (2009) have 
inferred an even higher value of OPR= 1.5-3.5 in the central 
molecular zone of the Galaxy. 

Using Herschel-HIFI observations of the fundamental rota- 
tional transitions of both ortho- and para-ammonia, it is for the 
first time possible to estimate the ammonia OPR in cold and dif- 
fuse interstellar gas of low excitation. Our results do, however, 
point to the surprising result of an OPR lower than unity. 

This is shown in Figs. 12 and 13 where the upper panels 
show the normalised lo-Oo ortho-NH3 and 2 1 — 1 1 para-NH3 
spectra, and the lower panels show the corresponding optical 
depth ratios for absorptions larger than 3cr as a function of LSR 
velocity towards both sources (channel widths are 0.26 km s _1 ). 
An ortho-to-para optical depth ratio of 4.2 corresponds to a col- 
umn density ratio of unity using the RADEX conversion factors 
in Table 2. The resulting column density ratios are given by the 
right hand y-axis in both figures. The para line has lower S/N 
than the ortho line in both sources and several velocity compo- 
nents are weak. Still, most column density ratios are found to be 
below unity within the uncertainties. 

We have also used Method I and II to estimate the OPR. 
Method I gives OPR(NH 3 ) = 0.7+0.1 and 0.6-0.7+0.1 in the 
strongest velocity components towards W49N (+39 km s _! com- 
ponent) and G10.6-0.4 (+16, 28, 30 and 41 km s" 1 compo- 
nents), respectively. The results from Method II gives similar 
results towards G10. 6-0.4: 0.5-0.7+0.1 in the velocity bins 
v LSR = 12.5-20, 25-29, 29-31 and 39.5-43.5 km s -1 ; but 
lower results towards W49N in v LSR = 36 - 42 km s -1 , 0.4+0.2. 

The errors are dominated by the noise and the uncertainty 
in the sideband gain ratio. We have used the errors in sideband 
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Fig. 13: G10.6-0.4 (W31C). Notation as in Fig. 12. 



gain ratios stated on the Herschel internet site (4% for ortho and 
6% for para). There is also an additional error arising from the 
rather large sideband separation of 12 GHz which means that the 
assumption of equal continuum temperature is not fully valid. 
This effect mimics a sideband gain ratio different from unity and 
is also taken into account (4 and 2% for ortho and para, respec- 
tively, estimated by 7c,l/7c,u = (vl/vuV 3 w ith /3 = 2). The er- 
rors from the noise is estimated by 5t = exp(r) x 51 Ik), and the 
errors from the sideband gain ratio and differences in the contin- 
uum in the sidebands are estimated by 5t — e - ln(l + eexp(r)) 
where e is both errors added in quadrature. The derivation of this 
formula assumes that the depth of the absorption line is correct, 
which means that 7c(SSB) - Ta is conserved but the continuum 
level is not. The noise and calibration errors are finally added in 
quadrature. The error estimates for Method I and II add respec- 
tive uncertainties from the methods to the calibration errors. 

Some previous measurements of ammonia have resulted 
in an OPR lower than unity, although not in diffuse gas. In 
Orion KL Hermsen et al. (1988) obtained OPR(NH 3 ) = 0.5 
from inversion emission lines, and Hermsen et al. (1985) found 
OPR( 15 NH 3 ) = 0.7. The latter low value was suggested to be 
caused by an excitation effect: a possible over-abundance of the 
unobserved K — state. This explanation cannot be applied in 
our case, since we actually observe the K = ortho state. Flower 
et al. (1995) modelled several NH 3 inversion transitions in the 
frame of C-shock models towards the warm and dense SgrB2 
envelope and their best OPR of ammonia was ~0.5, also lower 
than the statistical value of unity. Ceccarelli et al. (2002), how- 
ever, observed 2 1 high excitation ortho- and para-ammonia lines, 
both metastable and non-metastable levels, towards Sgr B2 using 
ISO, and derived an ammonia OPR of unity. 

There are to our knowledge no previous estimate of the 
OPR(NH 3 ) in diffuse gas, and it may very well be very differ- 
ent from the OPR in dense gas. We can, however, compare our 
strong absorption of the para-NH 3 line to the observations by 
Tieftrunk et al. (1994) who observed the two lowest inversion 
lines of para-NH 3 towards W49N. Their estimate of 7Y(p-NH 3 ) is 
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2.5xl0 12 crrT 2 which is less than half our value, ~6xl0 12 crrT 2 . 
If this column density is correct, it would imply an OPR higher 
than unity in this velocity component in W49N. In our excita- 
tion model, however, the integrated optical depth observed by 
Tieftrunk et al. (1994) in the 1 1 inversion line at +39 km s _1 im- 
plies a N(p-NH3) in good agreement with our value, since we 
obtain a higher excitation temperature than the 2.7 K used by 
Tieftrunk et al. (1994) for the 1 1 inversion transition. 

We have no clear explanation yet of our surprising result and 
we can only speculate about its origin. Either there are to us 
unknown instrumental effects, or that our assumption that the 
line-of-sight gas completely covers the background continuum 
within the beam is not correct, or there must exist some physi- 
cal or chemical processes that affect the OPR in diffuse gas. In 
general, the fact that interstellar clouds are weak plasmas allows 
for conversion processes at low temperatures that would not be 
present in a purely neutral gas. What OPR these processes lead 
to must be determined by careful modelling. In diffuse molecu- 
lar clouds of relatively low density, n(H2) < 10 3 , but relatively 
high ionisation fraction, n(e)/n(H2) ^ 10~ 4 , the rate of destruc- 
tion of NH3 by reactions with C + , H + , and , can be as high as 
10~ 9 s _1 at kinetic temperatures T ~ 30 K. The latter two ions can 
also interchange ortho and para states directly by proton substi- 
tution, perhaps at a comparable rate. These rates of destruction 
and interchange can approach the rates of radiative and colli- 
sional excitation out of the metastable levels, which means that 
the excitation and chemistry of NH3 should be treated together 
in a self-consistent fashion. Under these conditions, NH3 is a 
non-equilibrium system and a "spin temperature" derived from 
an ortho-to-para ratio is not expected to be meaningful. 

The NH| molecule could also affect the ammonia OPR in a 
similar manner. The fastest gas-phase reaction sequence of am- 
monia is 

N + NH + ^ H2 NHj ^ H2 NH3 ^ H2 NH4 NH 3 , NH 2 ,(2) 

where N + is formed by cosmic ray ionisation, or by reactions of 
He + with N2 or CN which are formed by neutral-neutral reac- 
tions. Whether or not ammonia is produced on grains or in the 
gas, its destruction in the gas by protonating ions such as H3 also 
leads to NITj, which can have a variety of nuclear spin states de- 
pending on the overall nuclear spin of Hi. The ammonia OPR 
could then simply reflect the spin states of NH^ . 

In addition, some laboratory evidence exists that dissociative 
recombination reactions have different rate coefficients depend- 
ing upon the nuclear spin configuration of the molecular ion. 

Observations of the ortho-to-para ratio of NH3 could thus 
provide a valuable insight into the competing processes of for- 
mation and destruction, radiative and collisional excitation, and 
reactive interchange processes. 

7. Summary 

Our spectrally resolved rotational transitions of NH, 0-NH2, 
ortho- and para-NH3 along the sight-lines towards the high- 
mass star-forming regions W49N and G10. 6-0.4 show remark- 
able similarities of line profiles and abundances. We find sim- 
ilar abundances of all three species and a co-existence in dif- 
fuse or translucent interstellar gas with a high molecular frac- 
tion. The mean abundance of ortho-NIT; abundance is ~2xl0~ 9 
towards both sources. The mean ratios of all three meth- 
ods of the nitrogen hydrides in all velocity components, are 
JV(NH)/JV(o-NH 3 ) = 5.9 and 3.5, and JV(o-NH 2 )/JV(o-NH 3 ) = 2.4 
and 2.0, towards G10.6-0.4 and W49N, respectively. This is in 



sharp contrast to previous observations of the nitrogen hydrides 
in dark clouds where the ammonia abundances are found to be 
~ 1 00 times higher than X(NH), and ~ 1 - 300 times higher than 
X(NH 2 ). NH and 0-NH2 are found to be linearly correlated with 
0-NH3 at least for iV(o-NH3) < 5xl0 12 cm -2 which corresponds 
to a few Ay. Upper limits of Af(NH + ) in both sources indicate a 
N(NH + )/N(NH) ratio of <2 - 14 %, with a mean of <6 %. 

Linear correlations are also found for all three nitrogen hy- 
drides with respect to CH, CN and HNC, although CH displays a 
more loose correlation than the latter two species. The nitrogen 
hydrides also largely follow the absorption pattern in Doppler 
velocity space of HCO + and water, a species also known to trace 
regions of a high molecular fraction. 

We have obtained a surprisingly low ortho-to-para ratio of 
ammonia, « 0.5-0.7+0.1, in the strongest velocity components, 
which is below the high-temperature limit of unity. No clear ex- 
planation has been found. More observations are needed of both 
the rotational transitions and the inversion lines with ground- 
based facilities, to be able to make firm conclusions about the 
ammonia OPR in diffuse gas. 

We will continue to investigate the absorption lines in the 
sight-lines towards the other six PRISMAS sources. This will 
allow an analysis of the nitrogen chemistry at various galac- 
tic distances from the Galactic Centre. We will also use new 
Open Time 1 (OT1) Herschel-HIFl data of higher excitation 
lines to analyse the hot core sources which will be compared and 
contrasted with the diffuse interstellar gas. The ortho-to-para ra- 
tio of NH3 will also be further investigated both in the sources 
and in the diffuse gas, in addition to the OPR of NH 2 , for which 
new OT1 data in four of the PRISMAS sources will be analysed 
and compared to the ammonia OPR. 
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Appendix A: Herschel observations 
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Table A.l: Herschel OBSID's of the observed transitions analysed in this paper. 
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Notes. {a) Three different frequency settings of the LO were performed, with approximately 15 km s -1 between each setting in order to determine 
the sideband origin of the signals. {b) Observed in the same setting as CH 2 at 946 GHz. ic) The ortho-NH 3 2 - l u transition at 1214.859 GHz is 
observed in the same setting as para-NH 3 2 t - li at 1215.246 GHz 
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Appendix B: Emission line contamination of o-NH 2 
towards W49N 



identified as the hfs blend in the lower half of the same spin- 
rotation doublet of NO 2 Yl l/2 J = 9.5 e -> 8.5 C at an average 
rest frequency of 952.145441 GHz. The PRISMAS observations 
have also found three additional NO lines in W49N shown in 
Fig. B.2, each one consisting of unresolved hfs components, 
while no emission of NO is found in G10. 6-0.4. 

Since the two NO lines seen in Fig. B.3 have almost equal 
line strength, and are also observed with the same instrument 
in the same band, we have used the observed 952.145 GHz 
transition as a template to remove the interfering NO line at 
952.464 GHz from the 0-NH2 absorption. In order to do this we 
use 

J 1 ^ Tq 

Tnorm = Tp > (B.l) 

JC + 'NO 

to calculate the normalised SSB intensity T noim in K, where Tq 
is the SSB continuum, and T^o is the intensity of the NO line 
(Ta - ?c,dsb)- The model NO line, shown in Fig. B.l, is then 
moved to the velocity of the emission line. The resulting absorp- 
tion line spectrum of 0-NH2 towards W49N with removed NO 
emission is shown in green in Fig. B.3. 
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Fig. B.l: W49N: Model of the 952.464GHz NO line to 
be removed from the 0-NH2 absorption lines. The model 
is using the NO 952.145 GHz emission line observed at 
vlsr = +7.7kms~ 1 and only shifted in velocity to a vlsr of 
+43.5 km s~' corresponding to a v L sr of +7.7kms~ 1 for that 
line. 



The 0-NH2 absorption towards W49N is contaminated by an 
emission line in the same sideband from the source, a more com- 
plicated situation than an emission line from the other sideband. 
In Fig. B.3 an emission line is clearly visible around 47 km s _1 in 
the 0-NH2 spectrum. Here, the intensities have been normalised 
to the continuum in single sideband as Ta/Tq-1 assuming a side- 
band gain ratio of unity where T& is the observed intensity and 
Tq is the SSB continuum as measured in line-free regions in 
the spectra. We identify the emission line as a blend of three 
unresolved hfs components of NO 2 ni/2 J = 9.5^ — > 8.5^, 
F = 10.5 - 9.5, 9.5 - 8.5, and 8.5 - 9.5, at 952.464201 GHz 
(weighted mean frequency, cf. Varberg et al. 1999). The emis- 
sion line appearing near 142 km s in this figure is consequently 
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Fig.B.2: W49N: Double sideband WBS spectra of four NO 
emission lines from the source itself. Each NO line consists of 3 
hfs components, with quatum numbers 2 Ili/2 J = 6.5 e — > 5.5 e 
(651 GHz), J = 6.5 f -> 5.5 f (652 GHz), J = 7.5 e -» 6.5 e 
(752 GHz), and J = 9.5 e -» 8.5 e (952 GHz). The emission line 
at-18kms-' blended with NO 651.433 GHz is SO ll n -lli . 
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Fig. B.3: NO emission from W49N in the line-of-sight 0-NH2 ab- 
sorption.: Normalised SSB spectrum of 0-NH2. The emission 
line at 142 km s" 1 is from NO at vo =952.145 GHz. The NO 
vo = 952.464GHz emission feature is seen at 47kms~ 1 in the 
0-NH2 absorption. The two NO lines have similar line strengths. 
The red line shows the SSB normalised 0-NH2 line without re- 
moval of the NO-line, and the green shows the spectra with re- 
moved NO-line. 
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Appendix C: Hyperfine Structure components Tab l e c -3: Hyperfine structure components of ortho-NH 2 1 1,1 - 

Oo,o- 



Table C.l: Hyperfine structure components of NH 3 2 
AT=1 <- 0,7 = 2 <- 1. 
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-18.90 


0.023 


974558.07 


9.82e-4 


-24.52 


0.035 


974564.78 


7.67e-4 


-26.58 


0.055 


974574.43 


8.15e-4 


-29.55 


0.088 


974583.03 


2.59e-4 


-32.20 


0.019 


974607.78 


1.21e-4 


-39.81 


0.013 



Notes. (a) The velocity offset from the strongest hfs component at 
974478.38 MHz. {b) The sum of the relative intensities of the 21 hfs 
components is 3.75. 



Table C.2: Hyperfine structure components of NH 3 S 
N=l <-0,7=0 <- 1. 



Frequency 




Av° 


Rel. Intensity 


(MHz) 


(s- 1 ) 


(kms- 1 ) 


AulXgu b 
A u /(main)xgj,(main) 


946380.79 


2.40e-3 


30.100375183108216 


0.369 


946380.79 


9.58e-4 


30.100375183108216 


0.295 


946419.79 


3.66e-4 


17.746741564471662 


0.056 


946419.98 


8.99e-4 


17.688682017264846 


0.276 


946475.82 


3.25e-3 





1.0 


946509.24 


1.93e-3 


-10.586919289910648 


0.297 


946509.24 


1.21e-3 


-10.586919289910648 


0.372 


946527.48 


1.81e-3 


-16.363416629849031 


0.278 


946527.56 


1.84e-4 


-16.389231431253044 


0.057 



Notes. ( "' The velocity offset from the strongest hfs component at 
946475.82 MHz. {b) The sum of the relative intensities of the 9 hfs com- 
ponents is 3.66. 



Frequency 
tlvlrlZJ 


(s ) 


Av" 
tKm s ) 


Rel. Intensity 

A, ( /(main)xg (( (main) 


952435.66 


4.86e-6 


44.91 


0.0002 


952446.99 


1.33e-5 


41.34 


0.0006 


952463.69 


3.07e-5 


36.09 


0.002 


952490.73 


3.5e-3 


27.58 


0.079 


952502.06 


1.28e-3 


24.01 


0.029 


952503.09 


1.79e-3 


23.69 


0.081 


952514.42 


3.65e-3 


20.12 


0.164 


952528.90 


2.03e-3 


15.56 


0.046 


952533.03 


2.23e-4 


14.27 


0.010 


952540.23 


6.53e-3 


12.00 


0.147 


952542.21 


7.02e-3 


11.37 


0.474 


952549.73 


2.63e-3 


9.01 


0.178 


952560.41 


3.32e-3 


5.65 


0.149 


952562.12 


6.42e-3 


5.11 


0.289 


952571.74 


7.55e-3 


2.08 


0.340 


952573.46 


3.88e-3 


1.54 


0.174 


952577.11 


8.47e-3 


0.39 


0.570 


952578.35 


l.lle-2 





1.0 


952600.46 


1.59e-3 


-6.96 


0.072 


952615.49 


3.58e-3 


-11.69 


0.081 


952626.82 


2.73e-3 


-15.25 


0.061 


952627.84 


2.67e-3 


-15.57 


0.120 


952628.25 


3.35e-3 


-15.70 


0.226 


952639.17 


1.40e-3 


-19.14 


0.063 


952653.66 


9.66e-3 


-23.70 


0.022 


952655.64 


7.43e-4 


-24.32 


0.050 


952659.49 


4.40e-4 


-25.54 


0.020 


952664.99 


1.58e-3 


-27.27 


0.036 


952686.88 


2.97e-4 


-34.16 


0.013 


952698.21 


7.06e-5 


-37.72 


0.003 



Notes. (a) The velocity offset from the strongest hfs component at 
952578.35 MHz. {b) The sum of the relative intensities of the 30 hfs 
components is 4.50. 

Table C.4: Hyperfine structure components of ortho-NH 3 1 -0q. 



Frequency" 


Ay* 


Rel. Intensity" "' 


(MHz) 


(kms- 1 ) 




572.4971 


-1.0473 


0.1999 


572.4984 


0.0 


1.0 


572.5002 


0.5237 


0.6 



Notes. (a) From Cazzoli et al. (2009). {b) The velocity offset from the 
strongest hfs component at 572.4984 MHz. {c) The sum of the relative 
intensities of the three hfs components is 1.8. 
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Table C.5: Hyperfine structure components of para-NH3 2\-\\. 



Frequency 
^ivinz ) 


A 

is ; 


LXV 
firm 

( Kill S ) 


K.el. Intensity 

A H ;(main)xg„(main) 


1215.2443 


1.02e-2 


0.31 


0.536 


1215.2449 


3.39e-3 


0.16 


0.179 


1215.2453 


5.65e-3 


0.06 


0.179 


1215.2456 


1.36e-2 





1.0 


1215.2459 


3.77e-4 


-0.09 


0.012 


1215.2468 


7.53e-3 


-0.32 


0.238 



Notes. (fl) Coudert & Roueff (2006, 2009). (h> The velocity offset from 
the strongest hfs component at 1215.2456 MHz. (c) The sum of the rela- 
tive intensities of the six hfs components is 2.14. 
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Appendix D: Method II results 
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Table D. 1 : Results from Method II" . 













W49N 










VLSR 




i \ \ u-LysrL2} 


i \ \ U-Ly jn 3 } 


(All } 


IN ri/U-lNIl 3 


U-1N n 2 / LJ-1N JA 3 


A ) 




Yin NH,Y 


firm c -1 ^ 


(cm ) 


(cm" 2 ) 


(cm 2 ) 


(cm ; 












36.0 - 42.0 


l.lel3 


8.3el2 


4.1el2 


6.3el3 


2.7 


2.0 


6.1e-9 


4.6e-9 


2.3e-9 


50.0 - 57.0 


1.4el2 


1.4el2 


1.2el2 


6.6el3 


1.2 


1.2 


7.4e-10 


7.4e-10 


6.4e-10 


57.0-61.0 


7.8el2 


3.5el2 


1.9el2 


7.3el3 


4.1 


1.8 


3.7e-9 


1.7e-9 


9.1e-10 


61.0 - 70.0 


1 1 1 

1.0el3 


4.8el2 


3.0el2 


9.1el3 


3.3 


1.6 


3.8e-9 


1 O,. i\ 

1.8e-9 


1.2e-9 


Total 


3.0el3 


1.8el3 


1.0el3 


2.9el4 


3.0 


1.8 


3.6e-9 


2.2e-9 


1.2e-9 


Mean 


1 .be.il 


4.3el2 


z.oelz 


7.3el3 


3.0 


1.8 


j.oe-y 


2.2e-y 


1.2e-y 


Median 


9.0el2 


a i„n 

4.zel2 


zoelz 


7.0el3 


3.6 


1.7 


4.5e-9 


2.1e-9 


1.2e-9 










G10.6-0.4 










Vlsr 


iV(NH) 6 


^(o-NHz) 6 


W(o-NH 3 ) 


N(CH) 


NH/o-NH 3 


o-NH 2 /o-NH 3 


X(NH) C 


X(o-NH 2 ) c 


X(o-NH 3 ) c 


(kms- 1 ) 


(cm- 2 ) 


(cm" 2 ) 


(cm" 2 ) 


(cm" 2 ) 












19^ 90(1 


z..yc 1 .3 


1 4pl ^ 


^ 4pl ? 

J.H-C 1Z 


1.2el4 


5.4 


z.u 


8 Q 
O.JC-7 


4 1p Q 


1 ftp 


20.0 - 25.0 


2.7el3 


1.2el3 


3.8el2 


9.3el3 


7.1 


3.2 


1.0e-8 


4.5e-9 


1.4e-9 


25.0 - 29.0 


3.0el3 


1.5el3 


8.4el2 


1.2el4 


3.6 


1.8 


8.8e-9 


4.4e-9 


2.5e-9 


29.0-31.0 


2.8el3 


1.2el3 


6.9el2 


5.4el3 


4.1 


1.7 


1.8e-8 


7.8e-9 


4.5e-9 


31.0 - 35.0 


1.4el3 


5.0el2 


2.8el2 


7.1el3 


5.0 


1.8 


6.9e-9 


2.5e-9 


1.4e-9 


35.0 - 39.5 


2.6el3 


8.2el2 


3.0el2 


l.lel4 


8.7 


2.7 


8.3e-9 


2.6e-9 


9.6e-10 


39.5 -43.5 


3.2el3 


l.lel3 


3.9el2 


6.3el3 


8.2 


2.8 


1.8e-8 


6.1e-9 


2.2e-9 


43.5 - 50.0 


1.8el2 


1.0el2 


8.0ell 


1.9el3 


2.3 


1.3 


3.3e-9 


1.8e-9 


1.5e-9 


Total 


1.9el4 


7.8el3 


3.5el3 


6.5el4 


5.4 


2.2 


1.0e-8 


4.2e-9 


1.9e-9 


Mean 


2.3el3 


9.8el2 


4.4el2 


8.1el3 


5.4 


2.2 


1.0e-8 


4.2e-9 


1.9e-9 


Median 


2.7el3 


l.lel3 


3.9el2 


8.2el3 


7.1 


3.0 


1.2e-8 


4.9e-9 


1.6e-9 



Notes. {a) Column densities, N(x), column density ratios, and abundances, X(x) = A'(x)/A'(H 2 ), in different velocity ranges. {b) Using ortho-NH 3 
1 -0 as a template. {c) Using [CH]/[H 2 ] = 3.5x 10~ 8 (Sheffer et al. 2008). 
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Appendix E: Method III results 
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Table E. 1 : Method III results" . 



W49N 



VLSR 


Av 


AT(NH) 


N(o-NH 2 ) 


jV(o-NH 3 ) 


A'(CH) 


NH/0-NH3 


0-NH2/0-NH3 


X(NH) 6 


X(o-NH 2 )* 


^(o-NHs) 6 


(kms- 1 ) 


(km s" 1 ) 




fcrrr 2 ) 


\ym } 


fcrrr 2 ) 












33.5 


1.9 


l.lel2 


2.2el2 


3.6ell 


3.2el3 


3.1 


6.1 


1.2e-9 


2.4e-9 


3.9e-10 


39.6 


1.1 


1.2el3 


7.4el2 


3.2el2 


3.3el3 


3.8 


2.3 


1.3e-8 


7.8e-9 


3.4e-9 


59.4 


3.0 


l.lel3 


3.5el2 


1.7el2 


7.0el3 


6.5 


2.1 


5.5e-9 


1.8e-9 


8.5e-10 


62.8 


2.4 


9.7el2 


3.2el2 


1.9el2 


5.8el3 


5.1 


1.7 


5.9e-9 


1.9e-9 


l.le-9 


Total: 




3.4el3 


1.6el3 


7.2el2 


1.9el4 


4.7 


2.3 


6.1e-9 


3.0e-9 


1.3e-9 


Mean: 




8.5el2 


4.1el2 


1.8el2 


4.8el3 


4.7 


2.3 


6.1e-9 


3.0e-9 


1.3e-9 


Median: 




1.0el3 


3.4el2 


1.8el2 


4.6el3 


5.8 


1.9 


8.0e-9 


2.6e-9 


1.4e-9 












run 6- 


-0 4 










^LSR 


Av 


Af(NH) 


2V(o-NH 2 ) 


/^(o-NHj) 


A'(CH) 


NH/0-NH3 


o-NH 2 /o-NH 3 


X(NH) fc 


X(o-NH 2 )* 


^(o-NHj)* 


(kms- 1 ) 


(kms- 1 ) 


(cm- 2 ) 


(cm- 2 ) 


(cm- 2 ) 


(cm- 2 ) 












tahlpc 1 ft ^ 
LclUlCS IU.j 


1 Q 


? 4p1 ^ 


Q ftpl 7 


^ dpi 9 


7 7p1 % 


7 1 


7 8 
z.o 


1 1p 8 


4 4p Q 


1 Sp Q 

1 .JC-7 


1Q 

17.U 


1 3 


8p1 7 


2.3el2 


U.UC1 1 


Z.UC1 J) 




3.0 


9 7p-9 


1 lp-9 
3. 1C-7 


8 1p-10 

0. 1C-1U 


22.2 


4.2 


1.9el3 


7.4el2 


1.7el2 


9.0el3 


11 


4.4 


7.4e-9 


2.9e-9 


6.6e-10 


22.9 


1.0 


l.lel3 


3.0el2 


1.0el2 


1.2el3 


11 


3.0 


3.2e-8 


8.8e-9 


2.9e-9 


Zj.U 


3.3 


C 1 7 


A Qa 1 7 

4.oeiz 


i. leiz 


4.uei 3 


7 ^ 

/.J 


A A 
4.4 


1 'Xa Q 


A 7a O 

4.ze-v 


q i^o 1 n 


ZO.U 


1 O 


o.zelj 


t.zet j 


o.4etz 


y.oel 3 


< n 
j.U 


1 O 


1 1 a Q 

i.ie-o 


A 7a O 


z.je-y 


30.0 


1.5 


3.1el3 


1.4el3 


6.8el2 


6.1el3 


4.6 


2.1 


1.8e-8 


8.0e-9 


3.9e-9 


32.2 


1.7 


1.2el3 


3.1el2 


8.9ell 


2.5el3 


13 


3.5 


1.7e-8 


4.3e-9 


1.2e-9 


36.2 


3.8 


2.2el3 


7.3el2 


2.5el2 


9.4el3 


8.8 


2.9 


8.2e-9 


2.7e-9 


9.3e-10 


39.1 


1.2 


1.4el3 


3.2el2 


5.5ell 


4.6el3 


25 


5.8 


l.le-8 


2.4e-9 


4.2e-10 


41.0 


1.7 


3.5el3 


l.lel3 


3.1el2 


8.5el3 


11 


3.5 


1.4e-8 


4.5e-9 


1.3e-9 


45.1 


1.3 


3.7el2 


1.6el2 


3.8ell 


1.4el3 


9.7 


4.2 


9.3e-9 


4.0e.9 


9.5e-10 


Total: 




2.2el4 


7.9el3 


2.8el3 


6.7el4 


7.7 


2.8 


l.le-8 


4.2e-9 


1.5e-9 


Mean: 




1.8el3 


6.6el2 


2.4el2 


5.6el3 


7.7 


2.8 


l.le-8 


4.2e-9 


1.5e-9 


Median: 




1.7el3 


6.1el2 


1.4el2 


5.4el3 


12 


4.3 


l.le-8 


4.0e-9 


9.2e-10 



Notes. (n) LSR velocity, vlsr, line width, Av, column densities, A'(x), column density ratios and abundances, X(x) = N(x)/N(H 2 ), of each velocity 
component. ib) X(x) = N(x)/N(H 2 ) using [CH]/[H 2 ] = 3.5xlQ- 8 (Sheffer et al. 2008). 
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Appendix F: Results for CN and HNC 
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Table F. 1 : Method I and III results for CN and HNC . 



W49N 








Method I 






Method III 




V LSR 


Av 


N(CN) 


AT(HNC) 


CN/0-NH3 


HNC/0-NH3 


iV(CN) 


A'(HNC) 


CN/0-NH3 


HNC/0-NH3 


(kms" 1 ) 


(kms" 1 ) 


(cm- 2 ) 


(cm- 2 ) 






(cm- 2 ) 


(cm" 2 ) 






33.5 


1.9 


8.1el2 


4.3ell 


18 


0.93 


7.5el2 


4.5ell 


21 


1.3 


39.6 


1.1 


6.6el3 


2.9el2 


15 


0.67 


6.3el3 


2.9el2 


20 


0.91 




^ n 
i.U 




L. lelz 


1 / 


i.U 




2.2el2 


19 


1 "3 
L.J 


62.8 


2.4 


2.8el3 


1.5el2 


13 


0.68 


2.6el3 


1.7el2 


14 


0.89 


Total: 




1.4el4 


6.9el2 


15 


0.76 


1.3el4 


7.3el2 


18 


1.0 


Mean: 




3.4el3 


1.7el2 


15 


0.76 


3.2el3 


1.8el2 


18 


1.0 


Median: 




3.2el3 


1.8el2 


15 


0.84 


3.0el3 


2.0el2 


16 


1.1 


G10.6-0.4 








Method I 






Method III 




VLSR 


Av 


iV(CN) 


Af(HNC) 


CN/0-NH3 


HNC/0-NH3 


N(CN) 


A'(HNC) 


CN/0-NH3 


HNC/0-NH3 


(kms- 1 ) 


(kms" 1 ) 


(cm- 2 ) 


(cm- 2 ) 






(cm- 2 ) 


(cm- 2 ) 






16.3 


1.9 


1.0el4 


3.6el2 


24 


0.86 


9.0el3 


4.1el2 


26 


1.2 


19.0 


1.3 


3.3el3 


1.3el2 


38 


1.5 


2.9el3 


1.5el2 


48 


2.5 


22.2 


4.2 


5.7el3 


2.7el2 


32 


1.5 


5.1el3 


3.1el2 


30 


1.8 


22.9 


1.0 


6.6el3 


2.5el2 


44 


1.7 


5.8el3 


2.9el2 


58 


2.9 


25.0 


3.3 


9.5el2 




6.8 




8.4el2 


1.7el2 


7.6 


1.5 


28.0 


1.9 


1.2el4 


4.1el2 


16 


0.56 


l.lel4 


4.7el2 


17 


0.73 


30.0 


1.5 


1.4el4 


5.8el2 


18 


0.73 


1.3el4 


6.7el2 


19 


1.0 


32.2 


1.7 


2.5el3 


1.0el2 


19 


0.77 


2.3el3 


1.2el2 


26 


1.3 


36.2 


3.8 


1.0el4 


5.1el2 


31 


1.6 


9.0el3 


5.9el2 


36 


2.4 


39.1 


1.2 


3.7el3 


l.lel2 


79 


2.3 


3.3el3 


1.3el2 


60 


2.4 


41.0 


1.7 


1.6el4 


5.9el2 


39 


1.4 


1.4el4 


6.8el2 


45 


2.2 


45.1 


1.3 


1.2el3 


8.5ell 


20 


1.4 


l.lel4 


9.7el2 


289 


26 


Total: 




8.6el4 


3.4el3 


25 


1.0 


8.7el4 


5.0el3 


31 


1.7 


Mean: 




7.2el3 


3.1el2 


25 


1.0 


7.4el3 


3.6el2 


53 


2.6 


Median: 




6.2el3 


2.7el2 


37 


1.5 


7.3el3 


4.1el2 


31 


1.7 



Notes. {a) LSR velocity, v L sr, line width, Av, and column densities, N, of each velocity component, and column density ratios. 



CM. Persson et al.: Nitrogen hydrides in interstellar gas II, Online Material p 13 
Table F.2: CN and HNC results from Method IF . 



W49N 



VLSR 


iV(CN)* 


yV(HNC)* 


CN/0-NH3 


HNC/o-NH 3 


(kms ') 


(cm ) 


(cm z ) 






JO.U - 42. U 




z.yeiz 


1 ^ 

LJ 


V. 1 1 


50.0 - 57.0 


4.4el2 


5.7ell 


3.7 


0.48 


57.0-61.0 


2.1el3 


1.7el2 


11 


0.89 


61.0-70.0 


2.6el3 


2.0el2 


8.7 


0.67 


Total 


1.0el4 


7.2el2 


10 


0.70 


Mean 


2.6el3 


1.8el2 


10 


0.70 


Median 


2.4el3 


1.9el2 


9.6 


0.76 






mo f\— n a. 






Vlsr 


iV(CN)* 


2V(HNC) 6 


CN/0-NH3 


HNC/0-NH3 


^KIll S ) 


(cm ) 


(cm ; 






12.5 - 20.0 


1.0el4 


5.8el2 


19 


1.1 


20.0 - 25.0 


9.1el3 


5.2el2 


24 


1.4 


25 - 29 


9.6el3 


5.6el2 


11 


0.67 


29.0-31.0 


8.8el3 


4.9el2 


13 


0.71 


31.0-35.0 


4.2el3 


2.3el2 


15 


0.82 


35.0 - 39.5 


8.0el3 


5.3el2 


27 


1.8 


39.5 -43.5 


1.4el4 


7.4el2 


36 


1.9 


43.5 - 50.0 


9.8el2 


l.lel2 


12 


1.4 


Total 


6.5el4 


3.8el3 


18 


1.1 


Mean 


8.1el3 


4.7el2 


18 


1.1 


Median 


9.0el3 


5.3el2 


23 


1.4 



Notes. {a) Column densities and ratios in different velocity ranges. {b) Using ortho-NH 3 1 -0 as a template. 
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Appendix G: Figures 




Velocity [km s ] 

Fig.G.l: G10.6-0.4 (W31C): Double sideband spectra of NH, 
ortho-NH 2 , ortho- and para- 14 NH 3 , ortho- 15 NH 3 , and NH + over 
the LSR velocity range -50 to 85 km s -1 . 
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Fig. G.2: Gl 0.6—0.4 (W31C): Normalised ammonia and the am- 
monia isotopologue 15 NH3. (The emission line at vlsr ~41 and 
72kms~ 1 is an SO2 line from the upper sideband.) 
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Fig. G.3: W49N: Fits of method II to NH and o-NH 2 , and their 
observed spectra. Residuals are found at the bottom. Template 
spectrum is 0-NH3 at 572 GHz. 
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Fig. G.4: W49N: Fits of method II to NH, o-NH 2 and 0-NH3, 
and their observed spectra. Residuals are found at the bottom. 
Template spectrum is CH at 532 GHz. 
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Method II, W49N: NH, o-NH 2 , 0-NH3 from CH 



2 - NH 974 GHz 




Residuals 

10 15 20 25 30 35 40 45 50 55 

Velocity [km/s] 



Fig.G.5: G10.6-0.4 (W31C): Fits of method II to NH and 
0-NH2, and their observed spectra. Residuals are found at the 
bottom. Template spectrum is 0-NH3 at 572 GHz. 
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Fig.G6: G10.6-0.4 (W31C): Fits of method II to NH, o-NH 2 
and 0-NH3, and their observed spectra. Residuals are found at 
the bottom. Template spectrum is CH at 532 GHz. 
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Fig. G.7: W49N: Method III (XCLASS) fits to the Nitrogen hy- 
drides and CH. 
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Fig.G.8: G10.6-0.4 (W31C): Method III (XCLASS) fits to the 
Nitrogen hydrides and CH. 
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Fig.G.9: Comparison of absorption lines towards W49N and G10.6-0.4. The intensities have been normalised to single sideband 
continuum. An NO emission line from the source (952.464GHz) has been removed from the 0-NH2 spectrum (see Sect. B.). 
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Fig. G.10: Comparison of absorption lines towards W49N and G10.6-0.4. The intensities have been normalised to single sideband 
continuum. 
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Appendix H: Comparison plots of different methods 

In Figs. H.l -H.4 we plot the resulting column densities from 
our three different methods. 

Note that since Method I and Method III use Gaussians with 
smaller line widths than the velocity bins of Method II towards 
G10. 6-0.4, the number of velocity components is not the same 
in this source. The +33 km s _1 component towards W49N is not 
used in Method II, which on the other hand models a velocity 
bin that is difficult to fit using Gaussians in Method I and III 
(50-57 km s _1 ). The comparison of the nitrogen hydrides vs. 
CH in W49N, in two of the four velocity bins, shows the only 
cases in which the results do not agree reasonably well between 
the methods. The resulting CH column densities also show a 
larger spread between the methods than the other species. The 
reason is not fully understood, but is probably caused by the 
different approaches to the CH modelling. Method II uses the 
deconvolved CH spectra as a template for the nitrogen hydrides 
and is thereby trying to fit the broader CH absorption to the more 
narrow features of the nitrogen hydrides with moderate success. 
The other two methods use an opposite approach: they use the 
output of the fitting of the nitrogen hydrides as an input to CH 
to fit Gaussians only in the same parts in velocity space, and are 
thereby trying to fit narrow features to the broader CH absorp- 
tion. The fit is good in some parts of the velocity space, and very 
bad in the velocity space in which CH has absorption but the 
nitrogen hydrides do not, which is expected. 
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Fig. H.l: Column density comparison plots. The results from all three methods are plotted. The scatter in the results can be consid- 
ered as an estimate of the errors. 
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Fig. H.2: Conf: Column density comparison plots 2. Notation as in Fig. H. 1 . 
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Fig. H.3: Cont: Column density comparison plots 3. Notation as in Fig. H. 1 . 
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Fig. H.4: Cont: Column density comparison plots 4. Notation as in Fig. H. 1 . 



